arXiv:hep-ph/0512377vl 30 Dec 2005 


TASI 2004 Lecture Notes on Higgs Boson Physics* 


Laura Reina^ 

^Physics Department, Florida State University, 

315 Keen Building, 

Tallahassee, FL 32306-4350, USA 
E-mail: reina@hep.fsu.edu 

Abstract 

In these lectures I briefly review the Higgs mechanism of spontaneous symmetry breaking and 
focus on the most relevant aspects of the phenomenology of the Standard Model and of the Minimal 
Supersymmetric Standard Model Higgs bosons at both hadron (Tevatron, Large Hadron Collider) 
and lepton (International Linear Collider) colliders. Some emphasis is put on the perturbative 
calculation of both Higgs boson branching ratios and production cross sections, including the most 
important radiative corrections. 


* These lectures are dedicated to Filippo, who listened to them before he was born and behaved really well 
while I was writing these proceedings. 
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I. INTRODUCTION 

The origin of the electroweak symmetry breaking is among the most important open 
questions of contemporary particle physics. Since it was hrst proposed in 1964 by Higgs, 



Kibble, Guralnik, Hagen, Englert, and Brout |l|,l2,|3|, the mechanism of spontaneous sym¬ 
metry breaking known as Higgs mechanism has become part of the Standard Model (SM) 
of particle physics and of one of its most thoroughly studied extensions, the Minimal Super- 
symmetric Standard Model (MSSM). A substantial theoretical and experimental effort has 
been devoted to the study of the physics of the single scalar Higgs boson predicted by the 


realization of the Higgs mechanism in the SM, and of the multiple scalar and pseudoscalar 
Higgs bosons arising in the MSSM. Indeed, the discovery of one or more Higgs bosons is 
among the most important goals of both the Tevatron and the Large Hadron Collider (LHC), 
while the precise determination of their physical properties strongly support the need for a 
future high energy International Linear Collider (ILC). 

In these lectures I would like to present a self contained introduction to the physics of 
the Higgs boson(s). Civen the huge amount of work that has been done in this held, I 
will not even come close to being exhaustively complete. This is not actually my aim. For 
this series of lectures, I would like to present the reader with some important background 
of informations that could prepare her or him to explore further topical issues in Higgs 
physics. Also, alternative theoretical approaches to the electroweak symmetry breaking and 
the generation of both boson and fermion masses will not be considered here. They have 
been covered in several other series of lectures at this school, and to them I refer. 

In Section HD after a brief glance at the essence of the Higgs mechanism, I will review 
how it is embedded in the Standard Model and what constraints are directly and indirectly 
imposed on the mass of the single Higgs boson that is predicted in this context. Among the 
extensions of the SM I will only consider the case of the MSSM, and in this context I will 
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mainly focus on those aspects that could be more relevant to distinguish the MSSM Higgs 
bosons. Section uni will review the phenomenology of both the SM and the MSSM Higgs 
bosons, at the Tevatron and the LHC, and will then focus on the role that a high energy 
ILC could play in this context. Finally, in Section HVl I will briefly summarize the state of 
the art of existing theoretical calculations for both decay rates and production cross sections 
of a Higgs boson. 

Let me conclude by pointing the reader to some selected references available in the litera¬ 
ture. The theoretical bases of the Higgs mechanism are nowadays a matter for textbooks in 
Quantum Field Theory. They are presented in depth in both Refs, jj] and j^. An excellent 
review of both SM and MSSM Higgs physics, containing a very comprehensive discussion 
of both theoretical andphenomenological aspects as well as an exhaustive bibliography, has 
recently appeared jy, l7|. The phenomenology of Higgs physics has also been thoroughly 


covered in a fairly recent review pap er 


proceedings and reports []T| 


previous summer schools 


UQ 


mMM 


w. 


lich can be complemented by several workshop 
KtI. h 71 h S| |. Finally, series of lectures given at 


can provide excellent references. 


II. THEORETICAL FRAMEWORK: THE HIGGS MECHANISM AND ITS CON¬ 
SEQUENCES. 


In Yang-Mills theories gauge invariance forbids to have an explicit mass term for the 
gauge vector bosons in the Lagrangian. If this is acceptable for theories like QED (Quantum 
Electrodynamics) and QCD (Quantum Chromodynamics), where both photons and gluons 
are massless, it is unacceptable for the gauge theory of weak interactions, since both the 
charged (W^) and neutral (Z°) gauge bosons have very heavy masses {Mw — 80 GeV, 
Mz — 91 GeV). A possible solution to this problem, inspired by similar phenomena happening 
in the study of spin systems, was proposed by several physicists in 1964 11,0,01, and it is 
known today simply as the Higgs mechanism. We will review the basic idea behind it in 
Section FlI Al In Section FlI HI we will recall how the Higgs mechanism is implemented in the 
Standard Model and we will discuss which kind of theoretical constraints are imposed on the 
Higgs boson, the only physical scalar particle predicted by the model. Finally, in Section Hi PI 
we will generalize our discussion to the case of the MSSM, and use its extended Higgs sector 
to illustrate how differently the Higgs mechanism can be implemented in extensions of the 
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SM. 


A. A brief introduction to the Higgs mechanism 

The essence of the Higgs mechanism can be very easily illustrated considering the case of 
a classical abelian Yang-Mills theory. In this case, it is realized by adding to the Yang-Mills 
Lagrangian 

= with = (aM" - , (1) 

a complex scalar field with Lagrangian 

= {DyyD.cP - V{(t>) = {DyyD.y - - A(0»' , (2) 

where and A > 0 for the scalar potential to be bounded from below. The 

full Lagrangian 

C = Ca + (3) 

is invariant under a U{1) gauge transformation acting on the helds as: 

0(x) —> ) A^{x) —> A^{x) -I- -d^a{x) , (4) 

9 

while a gauge held mass term (i.e., a term quadratic in the helds A^) would not be gauge 
invariant and cannot be added to C if the U{1) gauge symmetry has to be preserved. Indeed, 
the Lagrangian in Eq. can still describe the physics of a massive gauge boson, provided 
the potential V (0) develops a non trivial minimum (0*07^ 0). The occurrence of a non trivial 
minimum, or, better, of a non trivial degeneracy of minima only depends on the sign of the 
gf parameter in Y(0). For p^>0 there is a unique minimum at 0*0 = 0, while for /i^<0 the 
potential develops a degeneracy of minima satisfying the equation 0*0 = —gf/iflX). This 
is illustrated in Fig. [TJ where the potential Y(0) is plotted as a function of the real and 
imaginary parts of the held 0 = 0i -|- i02. In the case of a unique minimum at 0*0 = 0 the 
Lagrangian in Eq. © describes the physics of a massless vector boson (e.g. the photon, in 
electrodynamics, with g = —e) interacting with a massive charged scalar particle. On the 
other hand, something completely diherent takes place when gf <t). Choosing the ground 
state of the theory to be a particular 0 among the many satisfying the equation of the 
minimum, and expanding the potential in the vicinity of the chosen minimum, transforms the 
Lagrangian in such a way that the original gauge symmetry is now hidden or spontaneously 
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FIG. 1: The potential F((/>) (^= (/>i + i4>2) plotted for an arbitrary positive value of A and for an 
arbitrary positive (right) or negative (left) value of /r^. 


broken, and new interesting featnres emerge. To be more specific, let’s pick the following 
00 minimnm (along the direction of the real part of 0, as traditional) and shift the 0 held 
accordingly: 

/ 2 \ 1/2 . 

00 = ^ ^ ^ ^ 71 ^M^)) ■ (5) 

The Lagrangian in Eq. dSl) can then be rearranged as follows: 


C 



+ fiV; + l(d^it>2f + gvA„d'‘,t,2 + ... 

massive vector field massive scalar field Goldstone boson 


( 6 ) 


and now contains the correct terms to describe a massive vector held with mass m\ = g‘^v^ 
(originating from the kinetic term of C^), a massive real scalar held 0i with mass = —2p^, 
that will become a Higgs boson, and a massless scalar held 02, a so called Goldstone boson 
which conples to the gange vector boson A^. The terms omitted contain couplings between 
the 01 and 02 helds irrelevant to this discnssion. The gange symmetry of the theory allows 
ns to make the particle content more transparent. Indeed, if we parameterize the complex 
scalar held 0 as: 


= +H{x)) ^ -^{v + H{x)) , 


(7) 


the X degree of freedom can be rotated away, as indicated in Eq. dZD, by enforcing the 17(1) 
gange invariance of the original Lagrangian. With this gange choice, known as unitary gauge 
or unitarity gauge, the Lagrangian becomes: 


C = Ca + ^A^A^ + ^ [d^Hd^H + 2y^H^) + ... 


( 8 ) 
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which unambiguously describes the dynamics of a massive vector boson of mass m\ = 
g‘^v^, and a massive real scalar held of mass mjj = —2/i^ = 2Au^, the Higgs field. It is 
interesting to note that the total counting of degrees of freedom (d.o.f.) before the original 
U{1) symmetry is spontaneously broken and after the breaking has occurred is the same. 
Indeed, one goes from a theory with one massless vector held (two d.o.f.) and one complex 
scalar held (two d.o.f.) to a theory with one massive vector held (three d.o.f.) and one real 
scalar held (one d.o.f.), for a total of four d.o.f. in both cases. This is what is colorfully 
described by saying that each gauge boson has eaten up one scalar degree of freedom, 
becoming massive. 

We can now easily generalize the previous discussion to the case of a non-abelian Yang- 
Mills theory. Ca in Eq. (0 now becomes: 

= with F‘^^a,Al-a,Al + gr'-‘AlAl , (9) 

where the latin indices are group indices and are the structure constants of the Lie Al¬ 
gebra associated to the non abelian gauge symmetry Lie group, dehned by the commutation 
relations of the Lie Algebra generators = Let us also generalize the scalar 

Lagrangian to include several scalar helds fi which we will in full generality consider as real: 

- V (0) where V(f) = , (10) 

where the sum over the index i is understood and = df^ — igHA^. The Lagrangian of 
Eq. (j2I) is invariant under a non-abelian gauge transformation of the form: 

0i(x) ^ + , (11) 

A“(a:) ^ A“(a;) + ^d^a^{x) + /“''"Aj)(a;)a"(a:) . 

When /i^ < 0 the potential develops a degeneracy of minima described by the minimum 
condition: 0^ = (^q = —/r^/A, which only Exes the magnitude of the vector (pQ. By arbitrarily 
choosing the direction of fo, the degeneracy is removed. The Lagrangian can be expanded 
in a neighborhood of the chosen minimum and mass terms for the gauge vector bosons can 
be introduced as in the abelian case, i.e.: 

^ ... +tgHt'‘4,),{t'‘4,),A;A'’'‘ + ... ( 12 ) 

2 '-V-' 
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Upon diagonalization of the mass matrix in Eq. ()12|1 . all gauge vector bosons for 
which t°‘(f)Q 7 ^ 0 become massive, and to each of them corresponds a Goldstone particle, i.e. 
an unphysical massless particle like the y held of the abelian example. The remaining scalar 
degrees of freedom become massive, and correspond to the Higgs held H of the abelian 
example. 

The Higgs mechanism can be very elegantly generalized to the case of a quantum held 
theory when the theory is quantized via the path integral method^. In this context, the 
quantum analog of the potential U(0) is the effective potential Veff{(pci), dehned in term 
of the effective action r[0cz] (the generating functional of the IPI connected correlation 
functions) as: 

Veffi^Pci) = for (j)ci{x) = Constant = ipci , (13) 

where VT is the space-time extent of the functional integration and 4>ci{x) is the vacuum 
expectation value of the held conhguration ^(a;): 

0d(a;) = (H|0(x)|H) . (14) 


The stable quantum states of the theory are dehned by the variational condition: 




r[0, 


'cl\ 


= 0 


d 


'^effi^cl) — 0 , 


(15) 


<t>cl=‘Pcl difcl 

which identihes in particular the states of minimum energy of the theory, i.e. the stable 
vacuum states. A system with spontaneous symmetry breaking has several minima, all 
with the same energy. Specifying one of them, as in the classical case, breaks the original 
symmetry on the vacuum. The relation between the classical and quantum case is made 
even more transparent by the perturbative form of the ehective potential. Indeed, Veff{(pci) 
can be organized as a loop expansion and calculated systematically order by order in h: 


Veff{(pci) = V{ipci) -h loop ehects , (16) 

with the lowest order being the classical potential in Eq. ©• Quantum corrections to 
Veff{(pci) ahect some of the properties of the potential and therefore have to be taken into 


^ Here I assume some familiarity with path integral quantization and the properties of various generating 
functionals introduced in that context, as I did while giving these lectures. The detailed explanation of 
the formalism used would take us too far away from our main track 







account in more sophisticated studies of the Higgs mechanism for a spontaneously broken 
quantum gauge theory. We will see how this can be important in Section III (^1 when we 
discuss how the mass of the SM Higgs boson is related to the energy scale at which we 
expect new physics effect to become relevant in the SM. 

Finally, let us observe that at the quantum level the choice of gauge becomes a delicate 
issue. For example, in the unitarity gauge of Eq. m the particle content of the theory 
becomes transparent but the propagator of a massive vector held turns out to be: 




kP' — m\ 




m. 


(17) 


and has a problematic ultra-violet behavior, which makes more difficult to consistently dehne 
and calculate ultraviolet-stable scattering amplitudes and cross sections. Indeed, for the 
very purpose of studying the renormalizability of quantum held theories with spontaneous 
symmetry breaking, the so called renormalizable or renormalizability gauges gauges) are 
introduced. If we consider the abelian Yang-Mills theory of Eqs. (HD-©, the renormalizable 
gauge choice is implemented by quantizing with a gauge condition G of the form: 


1 

7!' 


G — -\- ^gv(j)2) , 


( 18 ) 


in the generating functional 


Z[J] = C J DAD(j)^D(j)2exp i J 


det 


Sa) 


(19) 


where C is an overall factor independent of the helds, ^ is an arbitrary parameter, and a 
is the gauge transformation parameter in Eq. ©• After having reduced the determinant in 
Eq. (fT!I|l to an integration over ghost helds (c and c), the gauge plus scalar helds Lagrangian 
looks like: 


C 


-G^ + C 


+ + ■ ■ 


+ C 


-d^ - iigv)^ 1 + 


(20) 


such that: 


(A^(fc)A^(-fc)) = 








k^k'' 


+ 




' k^k'' 
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( 21 ) 


—i 

{h{k)h{-k)) = {c{k)c{-k)) = 

where the vector held propagator has now a safe ultraviolet behavior. Moreover we notice 
that the 02 propagator has the same denominator of the longitudinal component of the 
gauge vector boson propagator. This shows in a more formal way the relation between the 
02 degree of freedom and the longitudinal component of the massive vector held A^, upon 
spontaneous symmetry breaking. 


B. The Higgs sector of the Standard Model 


The Standard Model is a spontaneously broken Yang-Mills theory based on the SU{2)l x 
U{1)y non-abelian symmetry groupj^ |5|. The Higgs mechanism is implemented in the 
Standard Model by introducing a complex scalar held 0, doublet of SU (2) with hypercharge 
= 1 / 2 , 


/ 


V 


( 22 ) 


with Lagrangian 

= (/^'^0)^/^^0 - /iVV - A(0V)' , (23) 

where D^cf) = — ig'Y^B^)^ and r“ = (t“/ 2 (for a = 1,2,3) are the SU{2) 

Lie Algebra generators, proportional to the Pauli matrix a°'. The gauge symmetry of the 
Lagrangian is broken to U{l)em when a particular vacuum expectation value is chosen, e.g.: 


( 0 ) = 


v/2 


with V = 


2\ 1/2 

^ ' (/i^ < 0, A > 0) 


A 


(24) 


Upon spontaneous symmetry breaking the kinetic term in Eq. (P?|) gives origin to the SM 
gauge boson mass terms. Indeed, specializing Eq. m to the present case, and using Eq. ( 121 , 
one gets: 

/ 


(DI‘4.)'D^4, ^ ... + -(0 ,.) {gAX + 9'B„) (g2l‘V'' + g'8“) 

In^ - 


+ ■ ■ ■ 


■ ■■ + 


2 4 L' 


9\Alf + g\Alf + {-gAl + g'B,y 


+ ■■■ 


(25) 
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( 26 ) 


One recognizes in Eq. (ESD the mass terms for the charged gauge bosons W^\ 




Mw = g-, 


(27) 


and for the neutral gauge boson 
^0 =_1 

while the orthogonal linear combination of and remains massless and corresponds to 
the photon held {A^)\ 

1 




v'F+S^ 


+ dBfj.) 


Ma = 0 , 


(28) 


the gauge boson of the residual 17(l)em gauge symmetry. 

The content of the scalar sector of the theory becomes more transparent if one works in 
the unitary gauge and eliminate the unphysical degrees of freedom using gauge invariance. 
In analogy to what we wrote for the abelian case in Eq. dZD, this amounts to parametrize 
and rotate the 0(a;) complex scalar held as follows: 


0(x) = 


“71“ 


/ 


0 


yV + H{x) 

after which the scalar potential in Eq. (ESI) becomes: 




/ 


0 


yV + H{x) 


(29) 


Cs = - XvH^ - -H^ = - y -ivih 

4 2^ \ 2 


MhH^ - -XH^ 


(30) 


Three degrees of freedom, the x“(a;) Goldstone bosons, have been reabsorbed into the lon¬ 
gitudinal components of the and weak gauge bosons. One real scalar held remains, 
the Higgs boson H, with mass M|f = —2/i^ = 2An^ and self-couplings: 


H, 


H' 


V- 


.. H = 


_3^Mi 

V 


H, 


H' 


N / 

N X 

X 

X s 
X s 


H 


^ H 




Furthermore, some of the terms that we omitted in Eq. (EH), the terms linear in the gauge 
bosons Wj^ and Z^, dehne the coupling of the SM Higgs boson to the weak gauge helds: 



1) «-X 



H 


= 


' H 
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We notice that the couplings of the Higgs boson to the gauge helds are proportional 
to their mass. Therefore H does not couple to the photon at tree level. It is important, 
however, to observe that couplings that are absent at tree level may be induced at higher 
order in the gauge couplings by loop corrections. Particularly relevant to the SM Higgs 
boson phenomenology that will be discussed in Section Iflfl are the couplings of the SM 
Higgs boson to pairs of photons, and to a photon and a weak boson: 



as well as the coupling to pairs of gluons, when the SM Lagrangian is extended through the 
QCD Lagrangian to include also the strong interactions: 



innnr s 


W(r(p 


The analytical expressions for the iLyy, H'yZ, and Hgg one-loop vertices are more in¬ 
volved and will be given in Section 1111 A1 As far as the Higgs boson tree level couplings 
go, we observe that they are all expressed in terms of just two parameters, either A and /i 
appearing in the scalar potential of £</, (see Eq. or, equivalently, Mh and v, the Higgs 
boson mass and the scalar held vacuum expectation value. Since v is measured in muon 
decay to be n = = 24:6 GeV, the physics of the SM Higgs boson is actually just 

function of its mass Mh- 

The Standard Model gauge symmetry also forbids explicit mass terms for the fermionic 
degrees of freedom of the Lagrangian. The fermion mass terms are then generated via gauge 
invariant renormalizable Yukawa couplings to the scalar held 0: 

V, - r«(5*y4 - r«L-y4 + h.c. (31) 

where cf)^ = , and Tj (/ = u,d,l) are matrices of couplings arbitrarily introduced to 

realize the Yukawa coupling between the held 0 and the fermionic helds of the SM. and 
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(where i = 1, 2, 3 is a generation index) represent quark and lepton left handed doublets 
of SU{2)l, while and are the corresponding right handed singlets. When the 

scalar helds 0 acquires a non zero vacuum expectation value through spontaneous symmetry 
breaking, each fermionic degree of freedom coupled to 0 develops a mass term with mass 
parameter 

r ^ 

»/ = . 

where the process of diagonalization from the current eigenstates in Eq. m to the corre¬ 
sponding mass eigenstates is understood, and E f are therefore the elements of the diagonal¬ 
ized Yukawa matrices corresponding to a given fermion /. The Yukawa couplings of the / 
fermion to the Higgs boson (?//) is proportional to T/: 


f 



As long as the origin of fermion masses is not better understood in some more general 
context beyond the Standard Model, the Yukawa couplings i/f represent free parameter of the 
SM Lagrangian. The mechanism through which fermion masses are generated in the Stan¬ 
dard Model, although related to the mechanism of spontaneous symmetry breaking, requires 
therefore further assumptions and involves a larger degree of arbitrariness as compared to 
the gauge boson sector of the theory. 



C. Theoretical constraints on the Standard Model Higgs boson mass 

Several issues arising in the scalar sector of the Standard Model link the mass of the 
Higgs boson to the energy scale where the validity of the Standard Model is expected to 
fail. Below that scale, the Standard Model is the extremely successful effective held theory 
that emerges from the electroweak precision tests of the last decades. Above that scale, the 
Standard Model has to be embedded into some more general theory that gives origin to a 
wealth of new physics phenomena. From this point of view, the Higgs sector of the Standard 
Model contains actually two parameters, the Higgs mass {Mh) and the scale of new physics 
(A). 
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In this Section we will review the most important theoretical constraints that are imposed 
on the mass of the Standard Model Higgs boson by the consistency of the theory up to a 
given energy scale A. In particular we will touch on issues of unitarity, triviality, vacuum 
stability, hne tuning and, hnally, electroweak precision measurements. 


1. Unitarity 


The scattering amplitudes for longitudinal gauge bosons (V^Vl —^ VlVl, where V = 
W^,Z^) grow as the square of the Higgs boson mass. This is easy to calculate using 
the electroweak equivalence theorem J, |J, valid in the high energy limit (i.e. for energies 
s = Q'^^My), according to which the scattering amplitudes for longitudinal gauge bosons 
can be expressed in terms of the scattering amplitudes for the corresponding Goldstone 
bosons, i.e.: 




(33) 


where we have indicated by a;® the Goldstone boson associated to the longitudinal component 
of the gauge boson H®. For instance, in the high energy limit, the scattering amplitude for 
satishes: 

'Ml 


A{WEWE WEWl) = A{uj^uj- + O 


i-W 


where 


A{ijJ^ 0 J —^ a;■''a; ) = 


M^ 


s-Ml 


+ 


t-Ml 


Using a partial wave decomposition, we can also write A as: 

CXD 

M = IbTT ^(2/ + 1)P;(cos 6*)a; , 


(34) 


(35) 


(36) 


1=0 


where a; is the spin I partial wave and Pi{cos6) are the Legendre polynomials. In terms of 
partial wave amplitudes ai, the scattering cross section corresponding to A can be calculated 
to be: 

IfiTT ^ 

(37) 


a = ^f;(2; + l)|a,p , 


1=0 


where we have used the orthogonality of the Legendre polynomials. Using the optical theo¬ 
rem, we can impose the unitarity constraint by writing that: 


°° 1 

<T =- Y.(2l + = -Im [A(e = 0)1 , 

^ 1=0 ^ 


(38) 
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where A{6 = 0) indicates the scattering amplitude in the forward direction. This implies 
that: 

\ai\^ = Re{aiY + Im(az)^ = Im(az) —> |Re(az)| < ^ . (39) 

Via Eq. (jSHI), different ai amplitudes can than provide constraints on Mh- As an example, 
let us consider the J = 0 partial wave amplitude oq for the —*• scattering we 

introduced above: 


ao — 


IbTTS J-s 


Adt = 




16?™^ 


2 + 




M\ 


H 


In the high energy limit {Mjj -C s), Oq reduces to: 


Oq 


Mjj<t:s Mjj 






(40) 


(41) 


from which, using Eq. dsni), one gets: 


Mh < 870 GeV . 


(42) 


Other more constraining relations can be obtained from different longitudinal gauge boson 
scattering amplitudes. For instance, considering the coupled channels like ZlZl, 

one can lower the bound to: 

Mh < 710 GeV . (43) 

Taking a different point of view, we can observe that if there is no Higgs boson, or equiv¬ 
alently if Mh 3> s, Eq. gives indications on the critical scale above which new 

physics should be expected. Indeed, considering again scattering, we see 

that: 


32^^' 

(44) 

from which, using Eq. we get: 


y/s^ < 1.8 TeV . 

(46) 


Using more constraining channels the bound can be reduced to: 


< 1.2 TeV . (46) 

This is very suggestive: it tells us that new physics ought to be found around 1-2 TeV, i.e. 
exactly in the range of energies that will be explored by the Tevatron and the Large Hadron 
Gollider. 
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2. Triviality and vacuum stability 


The argument of triviality in a theory goes as follows. The dependence of the quartic 
coupling A on the energy scale {Q) is regulated by the renormalization group equation 


dX{Q) 

dQ‘^ 



(47) 


This equation states that the quartic coupling A decreases for small energies and increases 
for large energies. Therefore, in the low energy regime the coupling vanishes and the theory 
becomes trivial, i.e. non-interactive. In the large energy regime, on the other hand, the 
theory becomes non-perturbative, since A grows, and it can remain perturbative only if A is 
set to zero, i.e. only if the theory is made trivial. 

The situation in the Standard Model is more complicated, since the running of A is 
governed by more interactions. Including the lowest orders in all the relevant couplings, we 
can write the equation for the running of X{Q) with the energy scale as follows: 

= 24A'^ - (3j'2 + 9/ - 24i;)X + + g/ - Hyf +■■■ (48) 


where t = lrL{Q‘^/Q^) is the logarithm of the ratio of the energy scale and some reference scale 
Qo square, yt = mt/v is the top-quark Yukawa coupling, and the dots indicate the presence 
of higher order terms that have been omitted. We see that when Mh becomes large, A 
also increases (since Mfj = 2Au^) and the hrst term in Eq. dominates. The evolution 
equation for A can then be easily solved and gives: 


A(Q) = 


X{Qz 


(49) 


1 - 4flA(<3(|) 111 (Ij) 

When the energy scale Q grows, the denominator in Eq. (Hi may vanish, in which case 
X{Q) hits a pole, becomes inhnite, and a triviality condition needs to be imposed. This is 
avoided imposing that the denominator in Eq. (jini) never vanishes, i.e. that X{Q) is always 
hnite and 1/X{Q) > 0. This condition gives an explicit upper bound on Mh'- 


Mh< 


31og(^) 


(50) 


obtained from Eq. inD by setting Q = A, the scale of new physics, and Qo = v, the electroweak 
scale. 
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On the other hand, for small Mh , i.e. for small A, the last term in Eq. (14811 dominates 
and the evolution of \{Q) looks like: 

A(A) = A(n)-^l/2log(^^^ . (51) 

To assure the stability of the vacuum state of the theory we need to require that A (A) > 0 
and this gives a lower bound for Mh'- 

TA /A2\ 

A(A)>0 — . (52) 

More accurate analyses include higher order quantum correction in the scalar potential and 
use a 2-loop renormalization group improved effective potential, 14 //, whose nature and 
meaning has been briefly sketched in Section III A1 


3. Indirect bounds from electroweak precision measurements 


Once a Higgs held is introduced in the Standard Model, its virtual excitations contribute 
to several physical observables, from the mass of the W boson, to various leptonic and 
hadronic asymmetries, to many other electroweak observables that are usually considered in 
precision tests of the Standard Model. Since the Higgs boson mass is the only parameter in 
the Standard Model that is not directly determined either theoretically or experimentally, it 
can be extracted indirectly from precision hts of all the measured electroweak observables, 
within the ht uncertainty. This is actually one of the most important results that can be 
obtained from precision tests of the Standard Model and greatly illustrates the predictivity 
of the Standard Model itself. All available studies can be found on the LEP Electroweak 
Working Group and on the LEP 


diggs Working Group Web pages [1^ |2i| as well as in their 


main publications |21|, |2^ 123, |2J, |2^ . An excellent recent series of lectures on the subject 


of Precision Electroweak Physics is also available from a previous TASI school j27|. 

The correlation between the Higgs boson mass Mh, the W boson mass Mw, the top- 
quark mass m^, and the precision data is illustrated in Figs. |21 and El Apart from the 
impressive agreement existing between the indirect determination of Mw and rut and their 
experimental measurements we see in Fig. |21 that the 68% GL contours from LEP, SLD, 
and Tevatron measurements select a SM Higgs boson mass region roughly below 200 GeV. 
Therefore, assuming no physics beyond the Standard Model at the weak scale, all available 
electroweak precision data are consistent with a light Higgs boson. 
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FIG. 2: Comparison of the indirect measurements of Mw and mt (LEP I+SLD data) (solid contour) 
and the direct measurement {pp colliders and LEP II data) (dashed contour). In both cases the 
68% CL contours are plotted. Also shown is the SM relationship for these masses as a function of 
the Higgs boson mass, mn- The arrow labeled Aa shows the variation of this relation if a{M‘^) is 


varied by one standard deviation. The left hand side plot is from Ref. m and corresponds to the 


Winter 2005 situation (see Eq. ESI)) the right hand side plot is from the Ref. 


2^ and corresponds 


to the Summer 2005 situation (see Eq. El). The comparison between Winter and Summer 2005 is 
shown in the right hand side plot. 


The actual value of Mh emerging from the electroweak precision fits strongly depends 
on theoretical predictions of physical observables that include different orders of strong and 
electroweak corrections. As an example, in Fig. El the magenta arrow shows how the yellow 
band would move for one standard deviation variation in the QED fine-structure constant 
a(mD. It also depends on the fit input parameters. As we see in Fig. El Mh grows for 
larger m*. The sensitivity of the indirect bound on Mh to is clearly visible both in Fig. El 
and in Fig. El where you can find the famous blue band plot. In both hgures, we compare 
the results of the Winter 2005 and Summer 2005 electroweak precision fits. As far as the 
Higgs boson mass goes, the main change between Winter and Summer 2005 has been the 


18 











FIG. 3: The 68% confidence level contour in mt and Mh for the fit to all data except the direct 
measurement of mj, indicated by the shaded horizontal band of Tier width. The vertical band 


shows the 95% CL exclusion limit on Mh from direct searches. From Ref. 


22 |. 


value of the top-quark mass. We go from: 

Mh = GeV 


Mh < 251 GeV (95% GL) 


for mt = 178. ±4.3 GeV . 


(63) 


in Winter 2005 


22^ . to: 


for mt = 172.7 ±2.9 GeV , 


(54) 


Mh = 91^^^ GeV 
Mh < 186 - 219 GeV (95% GL) 

in Summer 2005 0, We see in Fig. |21 that the overlap between the direct and indirect 
determination of Mh is greatly reduced when the value of m* decreases and at the same 
time the minimum of the band in Fig. |3] considerably shifts. While in the hrst case the 
electroweak precision hts are still largely compatible with the direct searches at LEP II that 
have placed a 95% GL lower bound on Mh at: 


Mh > 114.4 GeV , 


(55) 


in the second case a large region of the band in Fig. |3J in particular the region about 
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FIG. 4: Ax^ = ~ Xmin curve. The line is the result of the fit using all electroweak data; 

the band represents an estimate of the theoretical error due to missing higher order corrections. 
The vertical band shows the 95% CL exclusion limit on Mh from direct searches. The solid and 
dashed curves are derived using different evaluations of . The dotted curve includes 

low data. The left hand side plot is from Ref. 221 and corresponds to the Winter 2005 situation 


(see Eq. 1^ . the right hand side plot is from the Ref. 


2^ and corresponds to the Summer 2005 


situation (see Eq. El- Note the different horizontal scale in the two plots. 


the minimum, is already excluded, and values of Mh very close to the experimental lower 
bound seem to be favored. 

It is fair to conclude that the issue of constraining Mh from electroweak precision fits 
is open to controversies and, at a closer look, emerges as a not clear cut statement. With 
this respect. Fig. El illustrates the sensitivity of a few selected electroweak observables to the 
Higgs boson mass as well as the preferred range for the SM Higgs boson mass as determined 
from all electroweak observables . One can observe that Mw and the leptonic asymmetries 
prefer a lighter Higgs boson, while and the NuTeV determination of sin^ 9^ prefer a 
heavier Higgs boson. A certain tension is still present in the data. We could just think that 
things will progressively adjust and, after the discovery of a light Higgs boson at either the 
Tevatron or the LHC, this will result in yet another amazing success of the Standard Model. 
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FIG. 5: Preferred range for the SM Higgs boson mass Mh as determined from various electroweak 
observables. The shaded band shows the overall constraint on the mass of the Higgs boson as 


derived from the full data set. Prom Ref. 
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Or, one can interpret the situation depicted in Fig. El as an unavoidable indication of the 
presence of new physics beyond the Standard Model. Indeed, since the data compatible with 
a lighter Higgs boson are very solid, one could either interpret the data compatible with a 
larger value of Mh as an indication of new physics beyond the Standard Model, or one could 
drop them as wrong, and still, the Higgs boson mass would turn out to be so small not to 
be compatible anymore with the Standard Model, signaling once more the presence of new 
physics. 
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4- Fine-tuning 


One aspect of the Higgs sector of the Standard Model that is traditionally perceived as 
problematic is that higher order corrections to the Higgs boson mass parameter square con¬ 
tain quadratic ultraviolet divergences. This is expected in a theory and it does not pose 
a renormalizability problem, since a theory is renormalizable. However, although per se 
renormalizable, these quadratic divergences leave the inelegant feature that the Higgs boson 
renormalized mass square has to result from the adjusted or fine-tuned balance between 
a bare Higgs boson mass square and a counterterm that is proportional to the ultraviolet 
cutoff square. If the physical Higgs mass has to live at the electroweak scale, this can cause 
a hne-tuning of several orders of magnitude when the scale of new physics A (the ultraviolet 
cutoff of the Standard Model interpreted as an effective low energy theory) is well above 
the electroweak scale. Ultimately this is related to a symmetry principle, or better to the 
absence of a symmetry principle. Indeed, setting to zero the mass of the scalar helds in the 
Lagrangian of the Standard Model does not restore any symmetry to the model. Hence, the 
mass of the scalar helds are not protected against large corrections. 

Models of new physics beyond the Standard Model should address this hne-tuning prob¬ 
lem and propose a more satisfactory mechanism to obtain the mass of the Higgs particle(s) 
around the electroweak scale. Supersymmetric models, for instance, have the remarkable 
feature that fermionic and bosonic degrees of freedom conspire to cancel the Higgs mass 
quadratic loop divergence, when the symmetry is exact. Other non supersymmetric models, 
like little Higgs models, address the problem diherently, by interpreting the Higgs boson as 
a Goldstone boson of some global approximate symmetry. In both cases the Higgs mass 
turns out to be proportional to some small deviation from an exact symmetry principle, and 
therefore intrinsically small. 

As suggested in Ref. 


28| . the no fine-tuning condition in the Standard Model can be 
softened and translated into a maximum amount of allowed fine-tuning, that can be directly 
related to the scale of new physics. As derived in Section mu upon spontaneous breaking 
of the electroweak symmetry, the SM Higgs boson mass at tree level is given by Mfi = —2/i^, 
where is the coefficient of the quadratic term in the scalar potential. Higher order 
corrections to Mfi can therefore be calculated as loop corrections to i.e. by studying 
how the effective potential in Eq. (uni) and its minimum condition are modihed by loop 
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FIG. 6: The SM Higgs boson mass Mh as a function of the scale of new physics A, with all the 
constraints derived from unitarity, triviality, vacuum stability, electroweak precision fits, and the 


requirement of a limited fine-tuning. The em 
less than 1 part in 10 fine-tuning. From Ref. 


a region is consistent with all the constraints and 


corrections. If we interpret the Standard Model as the electroweak scale effective limit of a 
more general theory living at a high scale A, then the most general form of /i^ including all 
loop corrections is: 

CXD 

= y + A= y c„(A,) log"(A/Q) , (66) 

71=0 

where Q is the renormalization scale, Aj are a set of input parameters (couplings) and the Cn 
coefficients can be deduced from the calculation of the effective potential at each loop order. 
As noted originally by Veltman, there would be no hne-tuning problem if the coefficient of A^ 
in Eq. (jsni) were zero, i.e. if the loop corrections to /i^ had to vanish. This condition, known 
as Veltman condition^ is usually over constraining, since the number of independent c„ (set 
to zero by the Veltman condition) can be larger than the number of inputs A*. However 
the Veltman condition can be relaxed, by requiring that only the sum of a hnite number of 


23 









































































terms in the coefficient of is zero, i.e. requiring that: 

l^max 

^ c„(Ai)log"(A/MH) = 0 , (57) 

0 

where the renormalization scale /i has been arbitrarily set to Mh and the order n has been 
set to Umax, fixed by the required order of loop in the calculation of 14//. This is based on 
the fact that higher orders in n come from higher loop effects and are therefore suppressed 
by powers of (167r^)“^. Limiting n to nmax, Eq- can now have a solution. Indeed, if the 
scale of new physics A is not too far from the electroweak scale, then the Veltman condition 
in Eq. m can be softened even more by requiring that: 

'fT'max ^.2 

c„(A,) log”(A/AfH) < ^ . (58) 

This condition determines a value of Amax such that for A < Amax the stability of the 
electroweak scale does not require any dramatic cancellation in /i^. In other words, for 
A < Amax the renormalization of the SM Higgs boson mass does not require any fine-tuning. 
As an example, for Umax = 0, Cq = (327r^n^)“^3(2M^ -|- M| -|- Mjj — 4:mf), and the stability 
of the electroweak scale is assured up to A of the order of Anv ~ 2 TeV. For Umax = 1 the 
maximum A is pushed up to A ~ 15 TeV and for nmax = ‘^ up to A ~ 50 TeV. So, just going 
up to 2-loops assures us that we can consider the SM Higgs sector free of fine-tuning up to 
scales that are well beyond where we would hope to soon discover new physics. 

For each value of Umax, and for each corresponding Amax, Mh becomes a function of the 
cutoff A, and the amount of fine-tuning allowed in the theory limits the region in the (A, Mh) 
plane allowed to Mh{A). This is well represented in Fig.lHl where also the constraint from the 
conditions of unitarity (see Section ITl (1 1|1 . triviality (see Section ITl C 2|1 . vacuum stability (see 
Section Hi (1 2jl and electroweak precision fits (see Section Hi C 3jl are summarized. Finally, 
the main lesson we take away from this plot is that if a Higgs boson is discovered new physics 
is just around the corner and should manifest itself at the LHC. 

D. The Higgs sector of the Minimal Supersymmetric Standard Model 

In the supersymmetric extension of the Standard Model, the electroweak symmetry is 
spontaneously broken via the Higgs mechanism introducing two complex scalar SU{2)l 
doublets. The dynamics of the Higgs mechanism goes pretty much unchanged with respect 
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to the Standard Model case, although the form of the scalar potential is more complex and its 
minimization more involved. As a result, the and weak gauge bosons acquire masses 
that depend on the parameterization of the supersymmetric model at hand. At the same 
time, fermion masses are generated by coupling the two scalar doublets to the fermions via 
Yukawa interactions. A supersymmetric model is therefore a natural reference to compare 
the Standard Model to, since it is a theoretically sound extension of the Standard Model, 
still fundamentally based on the same electroweak symmetry breaking mechanism. 

Far from being a simple generalization of the SM Higgs sector, the scalar sector of a super- 
symmetric model can be theoretically more satisfactory because: (i) spontaneous symmetry 
breaking is radiatively induced {i.e. the sign of the quadratic term in the Higgs potential is 
driven from positive to negative) mainly by the evolution of the top-quark Yukawa coupling 
from the scale of supersymmetry-breaking to the electroweak scale, and (ii) higher order 
corrections to the Higgs mass do not contain quadratic divergences, since they cancel when 
the contribution of both scalars and their super-partners is considered (see Section IHCH) . 

At the same time, the fact of having a supersymmetric theory and two scalar doublets 
modihes the phenomenological properties of the supersymmetric physical scalar helds dra¬ 
matically. In this Section we will review only the most important properties of the Higgs 
sector of the MSSM, so that in Section uni we can compare the physics of the SM Higgs 
boson to that of the MSSM Higgs bosons. 

I will start by recalling some general properties of a Two Higgs Doublet Model in Sec¬ 
tion cind I will then specify the discussion to the case of the MSSM in Section FH D 21 

In Sections IHD^ and iHml T will review the form of the couplings of the MSSM Higgs 
bosons to the SM gauge bosons and fermions, including the impact of the most impor¬ 
tant supersymmetric higher order corrections. A thorough introduction to Supersymmetry 
and the Minimal Supersymmetric Standard Model has been given during this school by 
Prof. H. Haber to whose lectures I refer Q 


1. About Two Higgs Doublet Models 

The most popular and simplest extension of the Standard Model is obtained by consid¬ 
ering a scalar sector made of two instead of one complex scalar doublets. These models, 
dubbed Two Higgs Doublet Models (2HDM), have a richer spectrum of physical scalar helds. 
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Indeed, after spontaneous symmetry breaking, only three of the eight original scalar degrees 
of freedom (corresponding to two complex doublet) are reabsorbed in transforming the orig¬ 
inally massless vector bosons into massive ones. The remaining hve degrees of freedom 
correspond to physical degrees of freedom in the form of: two neutral scalar, one neutral 
pseudoscalar, and two charged scalar helds. 

At the same time, having multiple scalar doublets in the Yukawa Lagrangian (see Eq. (ED) 
allows for scalar flavor changing neutral current. Indeed, when generalized to the case of 
two scalar doublet 0^ and 0^, Eq. m becomes (quark case only): 


fc=l,2 


E rJ,tQl4>'=4 + iLc. , 

k=l,2 


(59) 


where each pair of fermions {i,j) couple to a linear combination of the scalar helds 0^ and 
0^. When, upon spontaneous symmetry breaking, the helds 0^ and 0^ acquire vacuum 
expectation values 

k 

^ for fc = 1,2 , (60) 

the reparameterization of Cmkawa of Eq. dSHD in the vicinity of the minimum of the scalar 
potential, with (for k = 1,2), gives: 

k k 

/^Yukawa = ^j,k^ “k “41] ^fj,k^ + ^-C. + FC Couplings , (61) 

k v2 v2 






M^. 

^3 


where the fermion mass matrices and are now proportional to a linear combination 
of the vacuum expectation values of 0^ and 0^. The diagonalization of M“- and Mfj does 
not imply the diagonalization of the couplings of the 0'^ helds to the fermions, and Flavor 
Changing (EC) couplings arise. This is perceived as a problem in view of the absence of 
experimental evidence to support neutral havor changing ehects. If present, these ehects 
have to be tiny in most processes involving in particular the hrst two generations of quarks, 
and a safer way to build a 2HDM is to forbid them all together at the Lagrangian level. This 
is traditionally done by requiring either that n-type and d-type quarks couple to the same 
doublet (Model I) or that w-type quarks couple to one scalar doublet while d-type quarks to 
the other (Model II). Indeed, these two diherent realization of a 2HDM can be justihed by 
enforcing on Cyukawa the following ad hoc discrete symmetry: 

f $1 ^ _$i and <1)2 ^ <1)2 


(62) 


-d* and > ±u^ 
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The case in which FC scalar neutral current are not forbidden (Model III) has also been 
studied in detail. In this case both up and down-type quarks can couple to both scalar 
doublets, and strict constraints have to be imposed on the FC scalar couplings in particular 
between the hrst two generations of quarks. 

2HDMs have indeed a very rich phenomenology that has been extensively studied. In 
these lectures, however, we will only compare the SM Higgs boson phenomenology to the 
phenomenology of the Higgs bosons of the MSSM, a particular kind of 2HDM that we will 
illustrate in the following Sections. 


2. The MSSM Higgs sector: introduction 


The Higgs sector of the MSSM is actually a Model H 2HDM. It contains two complex 
SU (2 )L scalar doublets: 



(63) 


with opposite hypercharge {Y = ±1), as needed to make the theory anomaly-free^. <Fi 
couples to the up-type and $2 to the down-type quarks respectively. Correspondingly, the 
Higgs part of the superpotential can be written as: 


Vh = (|/i|' + m2)|<Fi|2 + (|/i|2 + m2)|<F2|2-/xH6,,(<Fl«F^2 + h.c.) 


^2 I „/2 „ „2 

^ (i«Fii^-i$2r) +|-i$i$2r 


in which we can identify three different contributions 


(64) 


(i) the so called D terms, containing the quartic scalar interactions, which for the Higgs 
helds <Fi and $2 correspond to: 

(|$i|2 - |<F2p)' + y |<hl$2r , (65) 

with g and g' the gauge couplings of SU{2)l and U{1)y respectively; 


^ Another reason for the choice of a 2HDM is that in a supersymmetric model the superpotential should 
be expressed just in terms of superfields, not their conjugates. So, one needs to introduce two doublets 
to give mass to fermion fields of opposite weak isospin. The second doublet plays the role of in the 
Standard Model (see Eq. ED), where has opposite hypercharge and weak isospin with respect to 4>- 
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(ii) the so called F terms, corresponding to: 


ihr(i$ir+i$2r); 


( 66 ) 


(in) the soft SUSY-breaking scalar Higgs mass and bilinear terms, corresponding to: 


+ h.c.) . 


(67) 


Overall, the scalar potential in Eq. (El depends on three independent combinations of 
parameters, |/ip +m^, |/ip and jiB. One basic difference with respect to the SM case 

is that the qnartic coupling has been replaced by gauge couplings. This reduced arbitrariness 
will play an important role in the following. 

Upon spontaneous symmetry breaking, the neutral components of <hi and $2 acquire 
vacuum expectation values 



( 68 ) 


and the Higgs mechanism proceed as in the Standard Model except that now one starts 
with eight degrees of freedom, corresponding to the two complex doublets $1 and $ 2 - Three 


degrees of freedom are absorbed in making the and the massive. The W mass is 
chosen to be: = g^iyl + n|)/4 = and this hxes the normalization of vi and ^ 2 , 


leaving only two independent parameters to describe the entire MSSM Higgs sector. The 
remaining Eve degrees of freedom are physical and correspond to two neutral scalar fields 


= —(A/ 2 Re 02 ~ "^^ 2 ) sina + (A/ 2 Re 0 ° — Ui) coscK 
= {V2Rect)l-V2) cos a + (y 2 Re 0 ? - vi) sin a , 


(69) 


one neutral pseudoscalar held 



(70) 


and two charged scalar fields 


sin (3 + (j)f cos (3 , 


(71) 


where a and (3 are mixing angles, and tan/d = ni/n 2 - At tree level, the masses of the scalar 
and pseudoscalar degrees of freedom satisfy the following relations: 
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( 72 ) 




making it natural to pick Ma and tan (3 as the two independent parameters of the Higgs 
sector. 

Eq. (j72|l provides the famous tree level upper bound on the mass of one of the neutral 
scalar Higgs bosons, h^\ 

< M| cos 2/3 < M| , (73) 


which already contradicts the current experimental lower bound set by LEP H: Mh > 


93.0 GeV j26|. The contradiction is lifted by including higher order radiative corrections 
to the Higgs spectrum, in particular by calculating higher order corrections to the neutral 
scalar mass matrix. Over the past few years a huge effort has been dedicated to the calcula¬ 
tion of the full one-loop corrections and of several leading and sub-leading sets of two-loop 
corrections, including resummation of leading and sub-leading logarithms via appropriate 
renormalization group equation (RGE) methods. A detailed discussion of this topic can be 
found in some recent reviews a ISQI, IM{ and in the original literature referenced therein. For 
the purpose of these lectures, let us just observe that, qualitatively, the impact of radiative 
corrections on can be seen by just including the leading two-loop corrections propor¬ 

tional to y^, the square of the top-quark Yukawa coupling, and applying RGE techniques to 
resum the leading orders of logarithms. In this case, the upper bound on the light neutral 
scalar in Eq. dZSD is modihed as follows: 




log 




Y 

i2M|; 


(74) 


where M| = -|- Mf^)/2 is the average of the two top-squark masses, rrit is the running 

top-quark mass (to account for the leading two-loop QGD corrections), and Xt is the top- 
squark mixing parameter dehned by the top-squark mass matrix: 

/ 


^Qt + + Dl mtXt 

mtXt + ml + 


(75) 


with Xt = At—fi cot [3 {At being one of the top-squark soft SUSY breaking trilinear coupling), 
= (1/2 — 2/3 sin 6*vc)Af| cos 2/3, and = 2/3sin^ cos2/3. Fig. [7| illustrates the 

behavior of as a function of tan/3, in the case of minimal and maximal mixing. For large 
tan /3 a plateau (i.e. an upper bound) is clearly reached. The green bands represent the 
variation of Mh as a function of mt when mi = 175±5 GeV. If top-squark mixing is maximal. 
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FIG. 7: The mass of the light neutral scalar Higgs boson, h^, as a function of tan/3, in the 
minimal mixing and maximal mixing scenario. The green bands are obtained by varying the 
top-quark mass in the m* = 175 ± 5 GeV range. The plot is built by fixing Ma = 1 TeV and 
MgusY = Mq =Mu = M£) = 1 TeV. From Ref. 


the upper bound on is approximately ~ 135 GeV^. The behavior of both and 

Mh± as a function of Ma and tan/3 is summarized in Fig. |H1 always for the case of maximal 
mixing. It is interesting to notice that for all values of Ma and tan/3 the Mh>MJJ^°‘^. Also 
we observe that, in the limit of large tan/3, i) for Ma< — Ma and Mh — 

while a) for Ma> Mh — Ma and Mh — 


3. MSSM Higgs boson couplings to electroweak gauge bosons 

The Higgs boson couplings to the electroweak gauge bosons are obtained from the kinetic 
term of the scalar Lagrangian, in strict analogy to what we have explicitly seen in the case 
of the SM Higgs boson. Here, we would like to recall the form of the HiVV and HiHjV 
couplings (for Hi = and V = W^, Z°) that are most important in order to 

understand the main features of the MSSM plots that will be shown in Section HTTI 

First of all, the couplings of the neutral scalar Higgs bosons to both and can be 

^ This limit is obtained for mt = 175 GeV, and it can go up to ~ 144 GeV for mt = 178 GeV. 
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FIG. 8: The mass of the light (h^) and heavy {H^) neutral scalar Higgs bosons, and of the 
charged scalar Higgs boson {H^) as a function of the neutral pseudoscalar mass Ma, for two 
different values of tan/3 (tan/3 = 3,30). The top-quark mass is fixed to mt = 174.3 GeV and 
Msusy^Mq=Mu = Mj:i = 1 TeV. The maximal mixing scenario is chosen. From Ref. 


written as: 


ghvv = gvMv sin(/3 - , gnvv = gvMy cos(/3 - a)g^'' , (76) 


where gv = ‘2My/v, while the A^VV and H^VV couplings vanish because of CP-invariance. 
As in the SM case, since the photon is massless, there are no tree level 'j'jHi and 
couplings. 

Moreover, in the neutral Higgs sector, only the and H^A^Z^ couplings are allowed 

and given by: 


gcos{(3-a) gsm{(3-a) 

ghAz = —^^— [Ph-PAr ) gHAz = —^^— {.ph-paY ■> An 

2 cos Ow 2 cos Uw 

where all momenta are incoming. We also have several HiHjV couplings involving the charge 
Higgs boson, namely: 


Ph+h-z — ^^^‘^^w{Ph+— Ph-Y 1 

2 cos Uw 

gH+H-y = —ie{PH+— Ph-Y ) 
gHThw^ = T'i^cos{(3 — a){ph — PhY^ i 
gH^HW± = - a){pH - PhY^ > 

gHTAw± = ^{pa-PhY^ ■ 
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At this stage it is interesting to introduce the so called decoupling limit, i.e. the limit of 
Ma 3> Mz, and to analyze how masses and couplings behave in this particular limit. Mh± 
in Eq. (17^ is unchanged, while become: 

Mh ~ and Mh ^ Ml + Ml sin^ 2/? . (79) 


Moreover, as one can derive from the diagonalization of the neutral scalar Higgs boson mass 
matrix: 


cos^(/9 — a) = 


Ml{A4l-Ml) «J»Mi 




(80) 


Ml(M% - Ml) 

From the previous equations we then deduce that, in the decoupling limit, the only light 
Higgs boson is hP with mass M^ ~ while Mh — Mh± — Ma 3> Mz, and because 

cos(/9 — a) ^ 0 (sin(/? — a) —> 1)), the couplings of hP to the gauge bosons tend to the SM 
Higgs boson limit. This is to say that, in the decoupling limit, the light MSSM Higgs boson 
will be hardly distinguishable from the SM Higgs boson. 

Finally, we need to remember that the tree level couplings may be modihed by radiative 
corrections involving both loops of SM and MSSM particles, among which loops of third 
generation quarks and squarks dominate. The very same radiative corrections that modify 
the Higgs boson mass matrix, thereby changing the dehnition of the mass eigenstates, also 
affect the couplings of the corrected mass eigenstates to the gauge bosons. This can be 
reabsorbed into the dehnition of a renormalized mixing angle a or a radiatively corrected 
value for cos(/3 — a) (sin(/3 — a)). Using the notation of Ref. j^, the radiatively corrected 
cos(/9 — a) can be written as: 

'Ml sin A/3 


cos{/3 — a) = K 


2Ml 


+ 0 


'MT 

Mi. 


where 


K = 1 + 


5M^ 


11 


5M^ 


22 


5M^ 


12 


(81) 


(82) 


2 M|cos2/? sin 2/5 ’ 

and 5M.ij are the radiative corrections to the corresponding elements of the CP-even Higgs 
squared-mass matrix (see Ref. j^). It is interesting to notice that on top of the traditional 
decoupling limit introduced above {Ma S> Mz), there is now also the possibility that cos{f3— 
a) —0 if iC —^ 0, and this happens independently of the value of Ma- 
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MSSM Higgs boson couplings to fermions 


As anticipated, $1 and $2 have Ynkawa-type conplings to the np-type and down-type 
components of all SU{2)l fermion donblets. For example, the Ynkawa Lagrangian corre¬ 
sponding to the third generation of qnarks reads: 


C 


Y ukawa 


= -hi 




- hh 


bRiflbL - bR(j )2 tL + h.c 


(83) 


Upon spontaneons symmetry breaking Cyukawa provides both the corresponding qnark 


masses: 


, Vi V sin j3 , , V 2 , V cos /3 

m = ht-^ = ht—j^ and mt = hb-^ = hb ^ 

and the corresponding Higgs-qnark conplings: 

-yt = [sin(/9 — a) -t- cot [3 cos(/9 — a)] yt , 

/ 

-yb = [sin(/9 — a) — tan (3 cos{f3 — a)] yb 


Qhtt 

9hbb 


COS a 

sm [3 
sin a 


(84) 


(85) 


COS (3 
sin a 


gmt — 


9Hbb — 


sin (3 
cos a 


cos (3 
9Ati = cot f3yt 
9 


9H±tb — 


yt = [cos(/3 - a) - cot (3 sin{(3 - a)] yt , 
yb = [cos(/9 — a) + tan f3 sin(/? — a)] yb , 
9Abb = f P Ub 1 

[mtCot[3{l - 75 ) +mbtsin[3{l + 75 )] 


2 a/2 iff 

where yq = mq/v (for q = t,b) are the SM conplings. It is interesting to notice that in the 
Ma 3> Mz deconpling limit, as expected, all the conplings in Eq. (|H3) rednce to the SM 
limit, i.e. all H^, and conplings vanish, while the conplings of the light nentral Higgs 
boson, rednce to the corresponding SM Higgs boson conplings. 

The Higgs boson-fermion conplings are also modihed directly by one-loop radiative cor¬ 
rections (sqnarks-glnino loops for qnarks conplings and sl^ton-nentralino loops for lepton 
conplings). A detailed discnssion can be fonnd in Ref. j?, 0] and in the literatnre refer¬ 
enced therein. Of particnlar relevance are the corrections to the conplings of the third qnark 
generation. These can be parameterized at the Lagrangian level by writing the radiatively 
corrected effective Ynkawa Lagrangian as: 


reff 

ukawa 


= —e 




{hb + 6hb)bjiQf^2 + {bit + bht)tRQf^\ 


( 86 ) 


- +h.c. , 
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where we notice that radiative corrections induce a small coupling between $1 and down- 
type helds and between $2 and up-type helds. Moreover the tree level relation between hb, 
ht, mh and rrit are modihed as follows: 


hbV 


mb = —j= cos 1 -f - -h 


5hb Ahbtan/3\ _ hbV 


72 

htv 


mt = -^sin/9 ——h 


] ~ 72 

6ht Ahbtan/9\ _ htV 


= ^ cos^(l -h Aft) , 


(87) 




ht ht 

where the leading corrections are proportional to Ah;, and turn out to also be tan (3 enhanced. 
On the other hand, the couplings between Higgs mass eigenstates and third generation 
quarks given in Eq. dini) are corrected as follows: 

1 Aht 


cos a 

9htt ~: aht 


sin (3' 


1 + Ai ht 


(cot (3 -|- tana) 
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9hbb = “ oVb 
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9 Abb = tan Pyb 
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1 
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— Abj (1 — tan a cot P) 
(cot P -\- tan P) 


(1 -h Ab)sm P \ hb 
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^mt cot P 
1 + 


1 - 


1 Aht 

1 + At ht 


(cot P -\- tan P) 


(1 + 75 ) 


'a_a 


(1 - 75) 


( 88 ) 


(1 -h Ab) sin^ P hb 

where the last coupling is given in the approximation of small isospin breaking effects, since 
interactions of this kind have been neglected in the Lagrangian of Eq. 


III. PHENOMENOLOGY OF THE HIGGS BOSON 

A. Standard Model Higgs boson decay branching ratios 

In this Section we approach the physics of the SM Higgs boson by considering its branch¬ 
ing ratios for various decay modes. In Section FlI HI we have derived the SM Higgs couplings 
to gauge bosons and fermions. Therefore we know that, at the tree level, the SM Higgs 
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FIG. 9: SM Higgs decay branching ratios as a function of Mh- The blue curves represent tree-level 
decays into electroweak gauge bosons, the red curves tree level decays into quarks and leptons, the 
green curves one-loop decays. From Ref. P|. 



boson can decay into pairs of electroweak gauge bosons {H —> W~^W~, ZZ), and into pairs 
of quarks and leptons (II QQ, while at one-loop it can also decay into two photons 

(H — > 77 ), two gluons (H — > gg), or a ■yZ pair (H jZ). Fig. 121 represents all the decay 
branching ratios of the SM Higgs boson as functions of its mass M^j. The SM Higgs boson 
total width, sum of all the partial widths T(H XX), is represented in Fig. II 01 

Fig.inishows that a light Higgs boson (Mh < 130 — 140 GeV) behaves very differently from 
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a heavy Higgs boson {Mh > 130 — 140 GeV). Indeed, a light SM Higgs boson mainly decays 
into a bb pair, followed hierarchically by all other pairs of lighter fermions. Loop-induced 
decays also play a role in this region. H ^ gg is dominant among them, and it is actually 
larger than many tree level decays. Unfortunately, this decay mode is almost useless, in 
particular at hadron colliders, because of background limitations. Among radiative decays, 
H —>■ 'j'j is tiny, but it is actually phenomenologically very important because the two photon 
signal can be seen over large hadronic backgrounds. On the other hand, for larger Higgs 
masses, the decays to W~^W~ and ZZ dominates. All decays into fermions or loop-induced 
decays are suppressed, except H ^ ti for Higgs masses above the ti production threshold. 
There is an intermediate region, around ~ 160 GeV, i.e. below the W^W~ and ZZ 
threshold, where the decays into WW* and ZZ* (when one of the two gauge bosons is off- 
shell) become important. These are indeed three-body decays of the Higgs boson that start 
to dominate over the H ^ bb two-body decay mode when the largeness of the HWW or 
HZZ couplings compensate for their phase space suppression^. The different decay pattern 
of a light vs a heavy Higgs boson influences the role played, in each mass region, by different 
Higgs production processes at hadron and lepton colliders. 

The curves in Fig. El are obtained by including all available QGD and electroweak (EW) 
radiative corrections. Indeed, the problem of computing the relevant orders of QGD and 
EW corrections for Higgs decays has been thoroughly explored and the results are nowadays 
available in public codes like HDEGAY Q. which has been used to produce Fig. El Indeed 
it would be more accurate to represent each curve as a band, obtained by varying the 
parameters that enters both at tree level and in particular through loop corrections within 
their uncertainties. This is shown, for a light and intermediate mass Higgs boson, in Fig. ITTI 
where each band has been obtained including the uncertainty from the quark masses and 
from the strong coupling constant. 

In the following we will briefly review the various SM Higgs decay channels. Giving a 
schematic but complete list of all available radiative corrections goes beyond the purpose of 
these lectures. Therefore we will only discuss those aspects that can be useful as a general 
background. In particular I will comment on the general structure of radiative corrections 

^ Actually, even four-body decays, corresponding to i/ —> W*W*,Z*Z* may become important in the 
intermediate mass region and are indeed accounted for in Fig. El 
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FIG. 11: SM Higgs boson decay branching ratios in the low and intermediate Higgs boson mass 
range including the uncertainty from the quark masses mt = 178 ± 4.3 GeV, mf, = 4.88 ± 0.07 GeV, 
and mc = 1.64 ± 0.07 GeV, as well as from as(Vfz) = 0.1172 ± 0.002. From Ref. Q]. 

to Higgs decay and I will add more details on QCD corrections to H ^ QQ {Q = heavy 


quark). 



For a detailed review of QCD correction in Higgs decays we refer the reader to Ref. |33|. 
Ref. also contain an excellent summary of both QCD and EW radiative corrections to 
Higgs decays. 

1. General properties of radiative correetions to Higgs deeays 

All Higgs boson decay rates are modified by both EW and QCD radiative corrections. 
QCD corrections are particularly important for H QQ decays, where they mainly amount 
to a redefinition of the Yukawa coupling by shifting the mass parameter in it from the pole 
mass value to the running mass value, and for H —> gg. EW corrections can be further 
separated into: i) corrections due to fermion loops, ii) corrections due to the Higgs boson 
self-interaction, and in) other EW corrections. Both corrections of type (ii) and (Hi) are in 
general very small if not for large Higgs boson masses, i.e. for Mh 3> Mw- On the other 
hand, corrections of type (i) are very important over the entire Higgs mass range, and are 
particularly relevant for Mh 2mt, where the top-quark loop corrections play a leading 
role. Indeed, for Mh -C 2mt, the dominant corrections for both Higgs decays into fermion 
and gauge bosons come from the top-quark contribution to the renormalization of the Higgs 


37 








wave function and vacuum expectation value. 

Several higher order radiative corrections to Higgs decays have been calculated in the 
large rrit limit, specifically in the limit when -C 2mt. Results can then be derived 
applying some very powerful low energy theorems. The idea is that, for an on-shell Higgs 
field = the limit of small masses {Mh -C 2mt) is equivalent to a pu —^ 0 limit, in 
which case the Higgs couplings to the fermion fields can be simply obtained by substituting 

^ 1^1 + , ( 89 ) 

in the (bare) Yukawa Lagrangian, for each massive particle i. In Eq. fl89jl is a constant 
held and the upper zero indices indicate that all formal manipulations are done on bare 
quantities. This induces a simple relation between the bare matrix element for a process 
with {X Y H) and without {X —>■¥)& Higgs held, namely 

lim .4(X ^Y + H) = 1 ^ Y) . (90) 

ph^o ^ drrv( 

When the theory is renormalized, the only actual diherence is that the derivative operation 
in Eq. o needs to be modihed as follows 

_, ^ _T (91) 

• 8m» 1 + 7 .„. dm, 

where is the mass anomalous dimension of fermion /j. This accounts for the fact that 
the renormalized Higgs-fermion Yukawa coupling is determined through the Z 2 and 
counterterms, and not via the Hff vertex function at zero momentum transfer (as used in 
the Ph ^ 0 limit above). 

The theorem summarized by Eq. (TO|) is valid also when higher order radiative corrections 
are included. Therefore, outstanding applications of Eq. o include the determination of 
the one-loop Hgg and ifyy vertices from the gluon or photon self-energies, as well as the 
calculation of several orders of their QCD and EW radiative corrections. Indeed, in the 
mt —> cx) limit, the loop-induced Efyy and Hgg interactions can be seen as effective vertices 
derived from an effective Lagrangian of the form: 

+ 0(a.)) , (92) 

i/TT V 

where F^a)^y is the field strength tensor of QED (for the TTyy vertex) or QCD (for the 
Hgg vertex). The calculation of higher order corrections to the H —>■ and H ^ gg 
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decays is then reduced by one order of loops! Since these vertices start as one-loop effects, 
the calculation of the hrst order of corrections would already be a strenuous task, and any 
higher order effect would be a formidable challenge. Thanks to the low energy theorem 
results sketched above, QCD NNLO corrections have indeed been calculated. 


2. Higgs boson decays to gauge bosons: H —> W~^W , ZZ 


The tree level decay rate for H —> VV {V = Z) can be written as: 


T{H VV) = 


GpMfj 




1 — Tv + -Ty ) Pv ) 


(93) 


16\/27r 

where f3v = — ry, ry = AMy/Mjj, and 5w,z = ‘^, 1 - 

Below the W^W~ and ZZ threshold, the SM Higgs boson can still decay via three (or 
four) body decays mediated by WW* {W*W*) or ZZ* [Z*Z*) intermediate states. As we 
can see from Fig. El the off-shell decays H —> WW* and H ^ ZZ* are relevant in the 
intermediate mass region around Mh — 160 GeV, where they compete and overcome the 
H ^ bb decay mode. The decay rates for H —> VV* —> Vfifj {V = hF^, Z) are given by: 

' Mw' 


T{H^WW*) = 
T{H ZZ*) = 




H 

5127r^" \Mh 

..4 


(94) 
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^ 40 2 

7- 


2048(1 - 5^)%^ V' 3"'^' 9 \Mh. 

where svy = sin is the sine of the Weinberg angle and the function F{x) is given by 

( 47 IS 1 \ 

—-1—;r) —3fl — Qx‘^ + ln(a;) 

2 2 x^ J ^ ' 

l-8a:2 + 20a:^ (‘ix^-l\ 

3- ,. — arccos | —— | . (95) 
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160 


’w 


M. 


A/4a;2 — 1 


2 a:^ 


3. Higgs boson decays to fermions: H —>■ QQ, Hi 


The tree level decay rate for H ^ ff {f = Q,l, Q =quark, I =lepton) can be written as: 

r{H ^ //) = , (96) 

where (3f = ^1 — r/, r/ = Arrij/Mfj, and (W)*’^ = 1,3. QCD corrections dominate over 
other radiative corrections and they modify the rate as follows: 

3G M 

T{H QQ)qcd = [‘^QCD + Ai] , (97) 
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where At represents specifically QCD corrections involving a top-quark loop. Aqcd and At 
have been calculated up to three loops and are given by: 


Aqcz) = 1 + + (35.94 - 1.36A^f) i 
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TT 


TT 


(98) 


(164.14 - 25J7Nf + 0.26Af|.) 

2 




H) 


TT 
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Cis{M, 
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TT 


1.57 


2 , Mjj 
-In ^ 
3 rrit 


1 

2^9 Ml 




where agiMn) and mQ{MH) are the renormalized running QCD coupling and quark mass 


in the MS scheme. It is important to notice that using the MS running mass in the overall 
Yukawa coupling square of Eq. (j^ is very important in Higgs decays, since it reabsorbs 
most of the QCD corrections, including large logarithms of the form ln(M^/mQ). Indeed, 
for a generic scale /r, is given at leading order by: 


mQ{fi)LO = rngirng) ( ( 99 ) 

where bo and 70 are the first coefficients of the (3 and 7 functions of QCD, while at higher 
orders it reads: 

mQ(/i) = rfigimo) 1 (100) 

where, from renormalization group techniques, the function f{x) is of the form: 


f{x) = 1 


12 

25 


(101) 


1 [1 + 1.014a;+ .. 

12 

23 

.] for mc<li<mb , 

f{x) = 1 

77 

1 [1 + 1.175X + .. 

1 

7 

.] for mb<^i<mt , 


f{x) = 1 

(t) 

[1 + 1.398a; + ...] 

for fi>mt . 



As we can see from Eqs. (Unni) and (nnn), by using the MS running mass, leading and 
subleading logarithms up to the order of the calculation are actually resummed at all orders 
in as- 

The overall mass factor coming from the quark Yukawa coupling square is actually the 
only place where we want to employ a running mass. For quarks like the b quark this could 
indeed have a large impact, since, in going from fi ~ Mh to /r ~ rrib, rfinifi) varies by almost 
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a factor of two, making therefore almost a factor of four at the rate level. All other mass 
corrections, in the matrix element and phase space entering the calculation of the H —> QQ 
decay rate, can in hrst approximation be safely neglected. 


4- Loop induced Higgs boson decays: H ^^,^Z,gg 


As seen in Section nm the if 77 and H'^Z couplings are induced at one loop via both a 
fermion loop and a W-loop. At the lowest order the decay rate for ii —77 can be written 
as: 

Gpa^Mfj 


T{H ^ 77) = 


128^2: 


TT"' 




f 


( 102 ) 


where = 1,3 (for f = l^q respectively), Qj is the charge of the / fermion species, 
Tf = Arrij/Mjj, the function /(r) is dehned as: 

r > 1 


/M = 



• 2 1 

1 arcsm 


In iTT 

1 4 

l—yl—T 


r < 1 


(103) 


and the form factors A^ and A^ are given by: 

Af = 2r[l + (1 -r)/(r)] , 
^w(^) = - [2 + 3r + 3r(2 - r)/(r)] 

On the other hand, the decay rate for H —> ■yZ is given by: 


(104) 


r(ii ^ jz) = 




1 - 


Ml 


X! ^w) 

f 


(105) 


where r* = AMf/Mfj and A* = AMf/M^ {i = f,W), and the form factors AJ{t,X) and 
A^{t, A) are given by: 

(106) 


■' COS Ow 


A(p(r, A) = cos 6 


w 


1 H— ) tan^ 9 
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sG- 
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A (a A) 


+ 4 (^3 - tan^ 0^/) i2(A A)| , (107) 

where N/ and Qf are dehned after Eq. (uni, and // is the weak isospin of the / fermion 
species. Moreover: 

^i(aA) = ^^^J /(^) - /(A)] + X ^ ^^^2 [^(^) -f/(A)] , 


2(r-A) 2(r-A)2' 


(r - A)^ 
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2(t - A) 


[/(t) - /(A)l , 
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9{r) = 


a/t — 1 arcsin ^ 
Vt 

^ ^In i+v^ - 


1—\/l—T 


?7r 


r > 1 
r < 1 


(108) 


(109) 


while /(r) is defined in Eq. fflOOj) . QCD and EW corrections to both T{H —> 77 ) and 
r{H —A 7 Z) are pretty small and for their explicit expression we refer the interested reader 
to the literatnre |(tL l83l|. 

As far as H gg is concerned, this decay can only be indnced by a fermion loop, and 
therefore its rate, at the lowest order, can be written as: 


^ gg) = 


CpalMl 


q 


( 110 ) 


where Tq = /{t) is defined in Eq. fjl 08 j) and the form factor A^{t) is given in 

Eq. (I100|l . QCD corrections to H ^ gg have been calcnlated np to NNLO in the rrit 
00 limit, as explained in Section IIII A 11 At NLO the expression of the corrected rate is 
remarkably simple 


r(Ef ^ gg{g), ggg) = Tlo{H gg) 


1 + E{tq) 


a 


(Nl) 


TT 


where 


^ Ar|r<4m2 95 7 33 - 2Nf 

£(r«) — + —6— 


log 


M. 


(Ill) 


( 112 ) 


When compared with the fnlly massive NLO calcnlation (available in this case), the two 
calcnlations display an impressive 10% agreement, as illnstrated in Fig. El even in regions 
where the light Higgs approximation is not jnstified. This is actnally dne to the presence 
of large constant factors in the first order of QCD corrections. We also observe that the 
first order of QCD corrections has qnite a large impact on the lowest order cross section, 
amonnting to more than 50% of Tfo on average. This has been indeed the main reason to 
prompt for a NNLO QCD calcnlation of r{H —a gg). The resnlt, obtained in the heavy-top 
approximation, has shown that NNLO QCD corrections amonnt to only 20% of the NLO 
cross section, therefore pointing to a convergence of the r{H —a gg) pertnrbative series. We 
will refer to this discnssion when dealing with the gg ^ H prodnction mode, since its cross 
section can be easily related to T(if —a gg). 
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FIG. 12: The QCD correction factor for the partial width T{H gg) as 
boson mass, in the full massive case with rrit = 178 GeV (dotted line) and 
limit (solid line). The strong coupling constant is normalized at as{Mz) = 


a function of the Higgs 
in the heavy-top-quark 
0.118. Prom Ref. 


B. MSSM Higgs boson branching ratios 


The decay patterns of the MSSM Higgs bosons are many and diverse, depending on the 
specihc choice of supersymmetric parameters. In particular they depend on the choice of 
Ma and tan /3, which parameterize the MSSM Higgs sector, and they are clearly sensitive to 
the choice of other supersymmetric masses (gluino masses, squark masses, etc.) since this 
determines the possibility for the MSSM Higgs bosons to decay into pairs of supersymmetric 
particles and for the radiative induced decay channels gg,'yy,'yZ) to receive 

supersymmetric loop contributions. 

In order to be more specific, let us assume that all supersymmetric masses are large enough 
to prevent the decay of the MSSM Higgs bosons into pairs of supersymmetric particles (a 
good choice could be Mg = Mq = Mu = M^ = 1 TeV). Then, we only need to examine 
the decays into SM particles and compare with the decay patterns of a SM Higgs boson to 
identify any interesting difference. From the study of the MSSM Higgs boson couplings in 
Sections III D 31 and III D 4L we expect that: i) in the decoupling regime, when Ma ^ the 
properties of the neutral Higgs boson are very much the same as the SM Higgs boson; 
while away from the decoupling limit ii) the decay rates of and to electroweak gauge 
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FIG. 13: Branching ratios for the and MSSM Higgs bosons, for tan P = 3, 30. The range 
of Mh corresponds to Ma = 90 GeV — 1 TeV, in the MSSM scenario discussed in the text, with 
maximal top-squark mixing. The vertical line in the left hand side plots indicates the upper bound 
on Mfi, which, for the given scenario is = 115 GeV (tan/3 = 3) or M^°'^ = 125.9 GeV 

(tan/3 = 30). From Ref. 


bosons are suppressed with respect to the SM case, in particular for large Higgs masses 
in) the —>• VV (V = W^, Z^) decays are absent, iv) the decay rates of hP and 

to T^T~ and bb are enhanced for large tan/3, v) even for not too large values of tan/3, due 
to ii) above, the h^,H^ t^t~ and h^,H^ — > bb decay are large up to the tt threshold, 
when the decay ^ tt becomes dominant, vi) for the charged Higgs boson, the decay 
dominates over th below the tb threshold, and vice versa above it. 

As far as QCD and EW radiative corrections go, what we have seen in Sections 1111 A 2f 
nil A 41 for the SM case applies to the corresponding MSSM decays too. Moreover, the truly 
MSSM corrections discussed in Sections III D 31 and III D 41 need to be taken into account and 
are included in FigsUHland d 


44 














































200 300 500 700 1000 

(GeV) 



niH* (GeV) 



niH* (GeV) 


FIG. 14: Branching ratios for the and MSSM Higgs bosons, for tan /3 = 3,30. The range 
of Mfj± corresponds to Ma = 90 GeV — _l_TeV, in the MSSM scenario discussed in the text, with 
maximal top-squark mixing. From Ref. 


C. Direct bounds on both SM and MSSM Higgs bosons 


LEP2 has searched for a SM Higgs at center of mass energies between 189 and 209 GeV. 
In this regime, a SM Higgs boson is produced mainly through Higgs boson strahlung from 
Z gauge bosons, Z* ^ HZ, and to a lesser extent through WW and ZZ gauge 

boson fusion, —>• WW,ZZ HugUg, He~^e~ (see Fig. USD. Once produced, it decays 

mainly into bb pairs, and more rarely into r+r” pairs. The four LEP2 experiments have 
been looking for: i) a. four jet hnal state {H ^ bb, Z ^ qq), ii) a missing energy hnal state 
{H —>■ bb, Z ^ vV), in) a leptonic hnal state {H bb, Z ^ Hl~) and iv) a specihc r-lepton 
hnal state {H bh, Z —>■ plus H —>■ , Z qq). The absence of any statistical 

signihcant signal has set a 95% CL lower bound on the SM Higgs boson at 

Mhsm > 114.4 GeV . 

LEP2 has also looked for the light scalar (h°) and pseudoscalar (H°) MSSM neutral Higgs 
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FIG. 15: SM and MSSM neutral Higgs boson production channels at LEP2. 
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FIG. 16: 95% CL exclusion limits for MSSM Higgs parameters from LEP2: (M;i,tan/3) (left) and 
(MA,tan/3) (right). Both the maximal and no-mixing scenarios are illustrated, for M 5 = 1 TeV 
and mt = 179.3 GeV. The dashed lines indicate the boundaries that are excluded on the basis of 
Monte Carlo simulations in the absence of a signal. Erom Ref. [3|. 
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FIG. 17: Leading Higgs production processes at hadron colliders: gg H, qq ^ QqH, and 
qq WH, ZH. 


bosons. In the deconpling regime, when is very heavy and behaves like a SM Higgs 
bosons, only can be observed and the same bonnds established for the SM Higgs boson 
apply. The bonnd can however be lowered when niA is lighter. In that case, and can 
also be pair prodnced throngh e'''e“ —>■ Z —>■ h^A^ (see Fig. Hi). Combining the different 
production channels one can derive plots like those shown in Fig. where the excluded 
(M/j,tan/?) and(Myi, tan/?) regions of the MSSM parameter space are shown. The LEP2 
collaborations have been able to set the following bounds at 95% CL: 

Mh,A > 93.0 GeV , 

obtained in the limit when cos(/3 — a) ~ 1 (anti-decoupling regime) and for large tan /?. 
The plots in Fig. ^1 have been obtained in the maximal mixing scenario (explained in 
Section FllP 2j) . For no-mixing, the corresponding plots would exclude a much larger region 
of the MSSM parameter space. 


Finally, the LEP collaborations have looked for the production of the MSSM c 


larged 


3J. An 


Higgs boson in the associated production channel: e’''e ■j, Z* H^H 

absolute lower bound of 

Mh± > 79.3 GeV 

has been set by the ALEPH collaboration, and slightly lower values have been obtained by 
the other LEP collaborations. 


D. Higgs boson studies at the Tevatron and at the LHC 

The parton level processes through which a SM Higgs boson can be produced at hadron 
colliders are illustrated in Figs. ITTIandlTHl 

Figures El and 1201 summarize the cross sections for all these production modes as functions 
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FIG. 18: Higgs production with heavy quarks: sample of Feynman diagrams illustrating the two 
corresponding parton level processes qq,gg —>■ ttH, bbH. Analogous diagrams with the Higgs boson 
leg attached to the remaining top(bottom)-quark legs are understood. 


of the SM Higgs boson mass, at the Tevatron (center of mass energy: a/s = 1.96 TeV) and at 
the LHC (center of mass energy: a/s = 14 TeV). These hgures have been recently produced 
during the TeV4LHC workshop and contain all known orders of QCD corrections as 
well as the most up to date input parameters. We postpone further details about QCD 
corrections till Section HU while we comment here about some general phenomenological 
aspects of hadronic Higgs production. 

The leading production mode is gluon-gluon fusion, gg ^ H (see hrst diagram in Fig. ED. 
In spite of being a loop induced process, it is greatly enhanced by the top-quark loop. For 
light and intermediate mass Higgs bosons, however, the very large cross section of this process 
has to compete against a very large hadronic background, since the Higgs boson mainly 
decays to bb pairs, and there is no other non-hadronic probe that can help distinguishing 
this mode from the overall hadronic activity in the detector. To beat the background, one 
has often to employ subleading if not rare Higgs decay modes, like H —> 77 , and this dilutes 
the large cross section. For larger Higgs masses, above the ZZ threshold, on the other hand, 
gluon-gluon fusion together with H ^ ZZ produces a very distinctive signal, and make this 
mode a ^^gold-plated mode" for detection. For this reason, gg —>■ H plays a fundamental 
role at the LHC over the entire Higgs boson mass range, but is of very limited use at the 
Tevatron, where it can only be considered for Higgs boson masses very close to the upper 
reach of the machine {Mh — 200 GeV). 

Weak boson fusion {qq qqH, see second diagram in Fig. ED ^md the associated 
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SM Higgs production 



m. [GeV] 


FIG. 19; Cross sections 
1.96 TeV). From Ref. 


35|. 


or SM Higgs boson production processes at the Tevatron, Run II (y^ = 



[GeV] 


FIG. 20: Cross sections for SM Higgs boson production processes at the LHC {^/s = 14: TeV). From 

Ref. 3- 


production with weak gauge bosons {qq —>• WH,ZH, see third diagram in Fig. ITTjl have 
also fairly large cross sections, of different relative size at the Tevatron and at the LHC. 
qq WH,ZH is particularly important at the Tevatron, where only a relatively light 
Higgs boson {Mh < 200 GeV) will be accessible. In this mass region, gg ^ H,H —>77 
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is too small and qq qqH is suppressed (because the initial state is pp). On the other 
hand, qq qqH becomes instrumental at the LHC {pp initial state) for low and intermedi¬ 
ate mass SM Higgs bosons, where its characteristic hnal state configuration, with two very 
forward jets, has been shown to greatly help in disentangling this signal from the hadronic 
background, using different Higgs decay channels. 

Finally, the production of a SM Higgs boson with heavy quarks, in the two channels 
QQydd QQH (with Q = t,b, see Fig. fTK|) . is sub-leading at both the Tevatron and the 
LHC, but has a great physics potential. The associated production with ti pairs is probably 
too small to be seen at the Tevatron, given the expected luminosities, but will play a very 
important role for a light SM Higgs boson {Mh < 130 —140 GeV) at the LHC, where enough 
statistics will be available to fully exploit the spectacular signature of a ttH, H ^ bb hnal 
state. Moreover, at the LHC, the associated production of a Higgs boson with top quarks will 
offer a direct handle on the top-quark Yukawa coupling (see Section [111 D2fl . On the other 
hand, the production of a SM Higgs boson with bb pairs is tiny, since the SM bottom-quark 
Yukawa coupling is suppressed by the bottom-quark mass. Therefore, the bbH, H ^ bb 
channel is the ideal candidate to provide evidence of new physics, in particular of extension 
of the SM, like supersymmetric models, where the bottom-quark Yukawa coupling to one 
or more Higgs bosons is enhanced (e.g., by large 
kinematically well within the reach of the Tevatron 
and D0 have already translated the absence of a bbh^, H^, signal into an upper bound 
on the tan/3 parameter of the MSSM. Were a signal observed, bbH could actually provide 
the hrst piece of evidence for new physics from RUN H. 

1. Searching for a SM Higgs boson at the Tevatron and the LHC 

Discovering a Higgs boson during RUN H of the Tevatron is dehnitely among the most 
important goal of this collider. It will be challenging and mainly luminosity limited, but 
recent studies have conhrmed that RUN H can push the 95% CL exclusion limit much farer 
than LEP2 and also shoot for a Sa or 5a discovery, depending on the integrated luminosity 
accumulated. 

The plot in Fig. |2I] shows the integrated luminosity that was originally estimated to be 
necessary to reach the 95% CL exclusion limit, the 3a, and the 5a discovery levels. It is 


tan/3 in the MSSM). bbH production is 
, RUN H. First studies from both CDF 
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FIG. 21: Integrated luminosity required for each experiment at the Tevatron, Run II, to exclude a 


SM Higgs boson at 95% CL or to observe it at the Scr or 5 it level. Original study from Ref. 


Q- 


given for a SM Higgs boson mass up to 200 GeV, that is to be considered as the highest 
Higgs boson mass reachable by RUN H. The curves have been obtained mainly by using the 
associated production with weak gauge bosons, qq VH (V = W^, Z^), with H ^ bb and 
H —> W~^W~, over the entire Higgs boson mass range, and gg ^ H with H ZZ in the 


upper mass region. As discussed in the introduction to Section UlIDl this can be understood 
in terms of production cross sections (see Fig. cni) and decay branching ratios (see Fig. El) 
over the Mh = 115 — 200 GeV mass range. From Fig. we see that with, e.g., 10 fb“^ of 
integrated luminosity RUN H will be able to put a 95% GL exclusion limit on a SM Higgs 
boson of mass up to 180 GeV, while it could claim a 3cr discovery of a SM Higgs boson 
with mass up to 125 GeV. A 5cr discovery of a SM Higgs boson up to 130 GeV, i.e. in the 
region immediately above the LEP2 lower bound, seemed to require 30 fb“^ of integrated 


luminosity, well beyond what is currently expected for RL 


N H. 


More recently, a new sensitivity study has appeared [1^, where the low mass region only 
has been revisited and new luminosity curves have been drawn. Mainly using the qq ^ VH 
{y = LF^, Z°) production mode, it appears that new analyses techniques will allow to obtain 
better results with less integrated luminosity. A 3 ct discovery of a SM Higgs boson with mass 
up to 125 GeV will now require only about 5 fb“^, while 10 fb“^ could allow a 5 ct discovery 
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of a SM Higgs boson with mass up to about 120 GeV, right at the LEP2 lower bound limit. 

At the LHC, all production modes will play an important role, thanks to the higher statis¬ 
tics available. In particular, it is natural to distinguish between a light {Mh < 130 — 140 GeV) 
and heavy {Mh > 130 — 140 GeV) mass region, as becomes evident by simultaneously looking 
at both production cross sections (see Fig. EDI) and decay branching ratios (see Fig. El) over 
the entire 115 —1000 GeV SM Higgs boson mass range. In the region of Mh < 130 —140 GeV 
the SM Higgs boson at the LHC will be searched mainly in the following channels: 

W+W-, ZZ , (113) 

qq^qqH ^ 77, W+W-, ZZ, r+r” , 

QQ, 99 ttH, H hi, r+r“ , 


while above that region, i.e. for Mh > 130 — 140 GeV, the discovery modes will be: 

gg^H,H ^ W+W-, ZZ , (114) 

qq qqH , H ^-f-f, W+W-, ZZ , 
qq,gg ^ ttH , H ^ IV+IV . 


These have been the modes used by both ATLAS and CMS to provide us with the discovery 
reach illustrated in Figs. E21 and ESI The ATLAS plots give the signal signihcance for a 
total integrated luminosity of 100 fb“^ (upper plot) and of 30 fb“^ (lower plot). The high 
luminosity (upper) plot belongs to the original ATLAS technical design report I 3 , and the 
weak boson fusion channels had not been studied in detail at that time. The lower luminosity 


(lower) plot is taken from a more updated study [^, and the weak boson fusion channels 
have been included in the low mass region, up to about Mh — 200 GeV, where they play 
an instrumental role towards discovery. Other instrumental channels in the low mass region 
are the inclusive Higgs production with iL —> 77 and, below Mj:/ = 130 GeV, tiH production 
with H — bl. In the high mass region, the inclusive production with H ZZ,WW 
dominates, although CMS has found a substantial contribution coming from weak gauge 
boson fusion with H ^ ZZ, WW. 
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FIG. 22: Significance for the SM Higgs boson discovery in various detection channels as a function 
of Mh- The upper plot is for 100 fb“^ of data and with no qqH channel included, the lower plot 
for 30 fb“^ and with the qqH channel included over the mass range Mh < 200 GeV. Results are 


from the ATLAS collaboration 


Q 


53 







FIG. 23: Luminosity required to reach a 5a discovery signal in CMS, using various detection 
channels, as a function of Mh- From Ref. B- 

2. Studies of a SM Higgs boson 

If a Higgs boson signal is established, the LHC will have the capacity of measnring several 
of its properties at some level of accuracy. In particular, it will be able to measure its mass, 
width, and couplings. At the same time, the charge and color quantum numbers of the 
newly discovered particle will be established by detecting a single production-decay channel; 
while a precise determination of its spin and parity will probably require more statistics 
than available at the LHC and will have to wait for a high energy Linear Collider to be 
established (see Section irrm . 

The two plots in Fig. show the precision with which both the mass and width of a 
SM-like Higgs boson will be determined by ATLAS and CMS combining 300 fb“^ of data 
per experiment. We see that below Mj:/~400 GeV the Higgs mass can be determined with 
a precision of about 0.1%, through H — ZZ 4/, complemented by H, WH, ttH{H 77) 
and ttH{H —> bb) in the low mass region. Above Mh—4:00 GeV the accuracy deteriorates for 
the smaller statistics available, although precisions of the order of 1% can still be obtained. 
We also see that the Higgs width above Mh — 200 GeV will be entirely determined through 
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FIG. 24: Expected precision in the measurement of a SM-like Higgs boson mass (M u) a nd width 
(F//) combining 300 fb~^ of data per experiment from ATLAS and CMS. From Ref. |4ri |. 

H ^ ZZ ^ 41, while below Mh — 200 GeV it is too small to be resolved experimentally 
and can only be determined indirectly, as we will discuss in the following. 

Finally, many studies in recent years have pointed to the fact that the LHC, under minimal 
theoretical assumptions, will have the potential to measure several Higgs boson couplings 
with an accuracy in the 10-30% range. The proposed strategy consists of measuring the 
production-decay channels listed in Eqs. ()113j) and ()114|) for a light {Mh < 130 - 140 GeV) 
or heavy {Mh > 130 — 140 GeV) Higgs boson respectively, and combine them to extract 
individual partial widths or ratios of partial widths. Indeed, if a given production-decay 
channel is observed, one can write that the experimentally measured product of production 
cross section times decay branching ratio corresponds, in the narrow width approximation, 
to the following expression: 

af{H) r,r„ 


{ap{H)Bi{H dd)y^P = 


Hth 

p 


( 116 ) 


where Fp and F^ are the partial widths associated with the production and decay channels 
respectively, while F is the Higgs boson total width. The coefficient (T*^(iL)/Fp^ can be 
calculated, while Fp, F^, and F is what needs to be determined. To each production-decay 
channel one can therefore associate a measurable observable 


yip) _ 


F„F 


pA d 


F 


(116) 
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where p and d label the production and decay channels respectively. is obtained from 
the experimental measurement of {ap{H)'Qi{H dd)Y^^, normalized by the theoretically 
calculable coefficient A signal in the {p,d) channel will measure zjY\ and 

therefore the product of Higgs couplings Ppy'^, since Fp ~ and Td — y^- Combining many 
different {p, d) channels, a system of equations of the form of Eq. (j1 Ibjl is obtained. Ratios 
of partial widths Fj/Fj, and therefore ratios of Higgs couplings, can then be derived in a 
model independent way, e.g.: 


F. 



(117) 


^ ^ y, 

r, ~ Vg ’ 

while individual partial width can be obtained with the further assumptions that: i) the total 
width is the sum of all SM partial widths, i.e. there is no new physics or invisible width 
contributions, and ii) the yw and yz couplings are related by the SU{2)l weak isospin 
symmetry. This is required by the fact that, in qq qqH, the W^W~ —> H and ZZ H 
fusion processes cannot be distinguished experimentally. 

The accuracy with which both individual couplings and ratios of couplings can be ex¬ 
tracted is mainly determined by the experimental error on the Z^^ measurements and by 
the theoretical uncertainty on the prediction of a^{H) in Eq. (Illbj) . Fig. OKI shows the 
estimated relative accuracy with which various channel will be detected at the LHC, assum¬ 
ing 200 fb“^ of data available for most channels. For the purpose of illustration. Fig. OKI 
shows the accuracy with which some of the SM Higgs partial width as well as its total width 
could be determined at the LHC, when the technique described above is implemented. The 
upper plots do not include any theoretical systematic error, while the lower plots include a 
theoretical systematic error of: 20% for gg H, 5% for qq qqH, and 10% for pp tiH. 
At the same time, all plots assumes the SU{2) gauge induced relation between y^ and yz, 
while the right hand plots also assume a SM-like relation between yb and yr- The more 
the assumptions, the better the accuracy with which the considered couplings can be deter¬ 
mined, and the more model dependence is introduced in the coupling determination. More 
sophisticated analyses have appeared in recent studies |l^ Q]- Overall, we can however 
conclude that the LHC has a great potential of giving a hrst fairly precise indication of the 
nature of the couplings of a Higgs boson candidate, although under some (well justihed) 
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FIG. 25: Relative accuracies on the measurement of the cross section of a scalar SM-like Higgs boson 
in the production+decay channels listed in Eqs. (inni) and (tTTH) . All channels have been rescaled to 
a total integrated luminosity of 200 fb~^, except pp ttH, H W~^W~ and qq ^ WH,H bb 
for which 300 fb“^ have been used, and gg ^ H,H W^W~ that was studies with 30 fb“^. 


From Ref. 


42| 


model assumptions. In particular, in the specific case of the top-quark Yukawa coupling, 
the LHC will be for a long time the only machine to be able to measure it with enough 
precision, since the measurement of yt in e^e“ —> ttH at a y^=500 GeV Linear Collider is 
statistically very limited, as we will see in Section [111 El 


3. Searching for a MSSM Higgs boson at the Tevatron and the LHC 

Most of the characteristics of the MSSM Higgs couplings that determine the pattern of 
decays reviewed in Section iTHRI also affect the mechanism of production of the MSSM Higgs 
bosons. In particular: 

• for Ma 3> Mz, the so called decoupling limit, 
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FIG. 26; Relative accuracy on the individual rates Fj expected at the LHC, when the ratio Fft/F,- is 
fixed to its SM value (right plots) or not (left plots). The integrated luminosity used is as explained 
in the caption of Fig. 1251 The upper plots show the accuracies obtained without including any 
theoretical systematic error, while the lower plots include the systematic theoretical errors explained 
in the text. From Ref. l42ll. 


hP —> Hsmi while 

Ma — Mh and g(^A,H)bb > dusMbl ) 9hvv < 9hsmVV, 


• while for Ma < Mz and tan (3^1: 

—^ 9hvv — 9hsmVV, while 

—^ Ma ^ Mh and 9i^A,h)bb > dHsubl ■, 9hvv < 9hsmVV- 

If we assnme snpersymmetric particles to be heavy enongh that the decay of a Higgs boson 
into snpersymmetric particles as well as the prodnction of a Higgs boson throngh the decay 
of a snpersymmetric particle is forbidden, the available prodnction processes are the SM 
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FIG. 27: Associated production modes for MSSM Higgs bosons. 


ones plus the associate production modes illustrated in Fig. 123 as well as the production of 
a charged Higgs boson via the decay of a t/t quark. A summary of the neutral MSSM Higgs 
boson production rates at the Tevatron and at the LHC is given in Figs. |2H1 and Eg for two 
values of tan/3, in the maximal mixing scenario (see Section FTTP 2jl . For the Tevatron the 
mass range is limited to the range kinematically accessible while for the LHC the entire Higgs 
mass range up to 1 TeV is covered. It is important to notice how for large tan (3 (tan /3 = 30 
in the plots of Figs. EHl and EH) the production of both scalar and pseudoscalar neutral 
Higgs boson with bottom quarks becomes dominant. In particular the inclusive production 
(denoted in the Figs. EHl and EHl as bb —>■ 0°, for 0° = H^, A°) becomes larger than the 

otherwise dominant gluon-gluon fusion mode {gg —>■ H), while the exclusive production 
(denoted in Figs. EHl and EHl as pp,pp —> bbcj)^) is right below gluon-gluon fusion but above 
all other production modes. More details on exclusive vs inclusive production of a Higgs 
boson with bottom quarks will be given in Section IIVI We also notice the subleading role 
played by vector boson fusion production {qq —> qqcf)^) due to the suppression (or absence in 
the case of A°) of the (jPVV couplings (for V = Z). Finally, the associated production 

modes for neutral Higgs bosons (see Fig. E3 have in general very small cross sections and 
are not considered in the plots of Figs. EHl and M 

For the case of the charged MSSM Higgs boson, we need to distinguish two cases: i) 
Mh± < rrit — mb and ii) Mh± > mt — m^- When Mh± < mt — mb, is mainly produced 
via the decay of a top or anti-top in tt production, i.e. via pp,pp tt with t bH^ or 
t bH~. On the other hand, when Mh± > mt — mb, is mainly produced through 
pp —>■ ibH^ ,tbH~, as well as through the tree level production modes qq —> H^H~ and 
bb , and at one loop through the associated modes gg H^H~ and gg —> . 

The overall cross section for the Tevatron and the LHC is illustrated in Figs. EDI and ED as 
a function of the charged Higgs boson mass, for different values of tan /3. The threshold 
behavior at mt — mb is clearly visible. 
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FIG. 28: Neutral MSSM Higgs boson cross sections at the Tevatron {^/s = 2 TeV). The upper 
plots are for the neutral scalar Higgs bosons, and while the lower plots are for the neutral 
pseudoscalar Higgs boson,The left hand plots are for tan/3 = 3 while the right hand plots are 
for tan/3 = 30. Prom Ref. 


The evidence or the absence of evidence for an MSSM Higgs boson at the Tevatron or 
at the LHC will place definite bounds on the parameter space of the MSSM. The reach 
of the Tevatron in the (Myi,tan/9) plane is illustrated in Fig. 1221 for different integrated 
luminosities. The shaded regions are spanned by using only the qq —>■ V(j)^ (0° —>■ bb) 
channel (with 0° = H^), while the region above the solid curves are covered by using the 

gg,qq —> 6&0° (0° —> bb) channel (with 0*^ = h^,H^,A^). The region below the black solid 
line is excluded by searches at LEP2. Fig. 1321 shows that, although discovery may require 
integrated luminosities that are beyond the reach of RUN2 of the Tevatron, with 5 fb“^ CDF 
and D0 will be able to exclude (in the maximal mixing scenario) almost all the parameter 
space of the MSSM at 95% C.L., a pretty impressive result by itself! Both CDF and D0 
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FIG. 29: Neutral MSSM Higgs boson cross sections at the LHC (\/s = 14 TeV). The upper hand 
plots are for the neutral scalar Higgs bosons, and while the lower plots are for the neutral 
pseudoscalar Higgs boson,^*^. The left hand plots are for tan/3 = 3 while the right hand plots are 
for tan/3 = 30. Prom Ref. 
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FIG. 31: Charged Higgs boson cross sections at the LHC as functions of Mfj±, for different values 
of tan p. From Ref. 


have indeed already presented results from searches conducted in the vv —^ bbcj)^ (0° — bb) 
channel with three or four 6-quark jets tagged in the hnal state 39|. The most recent 
results, from D0, are illustrated in Fig. 1221 where we can see that, depending on M^, values 
of tan f3 as low as tan /5 = 50 have already been excluded. The LHC 5cr discovery reach in 
the (M^,tan/3) parameter space is illustrated in Fig. Thanks to the high luminosity 
available and to the complementarity of various production and decay modes, the entire 
(M^,tan/?) parameter space can be covered, up to Ma of the order of 1 TeV. This gives us 
the exciting perspective that the LHC will be able to either discover or completely rule out 
the existence of an MSSM Higgs boson! 


E. Higgs boson studies at a future e"''e Linear Collider 

As we all know, an collider provides a very clean environment, with relatively 

simple signatures and very favorable signal to background ratios. Therefore, one of the most 
important roles that a high energy e’''e“ collider (to which we will refer as International 
Linear Collider (ILC) or simply Linear Collider (LC)) will play is to unambiguously identify 
any new particle discovered at the Tevatron or at the LHC, through a thorough program of 
precision measurements. This is true in Higgs boson physics as well, where we expect the 
mass, width, spin, and couplings of any Higgs boson candidates to be determined at the few 
percent level. In this Section I will mainly focus on a SM Higgs boson, since this is enough 
to illustrate the role played by a LC in Higgs physics, and since the true impact on the study 
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FIG. 32: 95% C.L. exclusion region and 5cr discovery region in the {Ma, tanP) plane, for the 
maximal mixing scenario. The shaded regions correspond to the search channel qq ^ V(fP {cjp 
bb) (with (j)^ = while the regions delimited by the colored solid lines correspond to the 

search channel gg,qq —> bhcfp {(jp bb) (with = hP,H^,A^). The two sets of lines correspond 
to simulations from CDF and D0 respectively. The region below the black line has been excluded 
by searches at LEP2. From Ref. See also Ref. 



FIG. 33: The 95% C.L. upper limit on tan/3 as a function of Ma for two scenarios of the MSSM, 
no mixing and maximal mixing. Also shown are the limits obtained by the LEP experiments for 


the same two scenarios of the MSSM. Erom Ref. 39|, D0 analysis based on 260 pb ^ of data. 
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FIG. 34: 5a discovery contours for MSSM Higgs boson detection in various channels, shown in 
the (M^,tan/3) parameter space, assuming the maximal mixing scenario. The left hand plot is 
from the CMS experiment and assumes 30 fb”^ of integrated luminosity, while the right hand plot 
combines ATLAS and CMS with a total integrated luminosity of 300 fb“^. 



FIG. 35: Higgs boson production via Higgs strahlung and W~^W , ZZ fusion at e"’“e colliders. 


of the MSSM parameter space will be shaped by the discoveries occured by the time a LC 

is built. 

First of all, the most important SM Higgs boson production processes in collisions 
are illustrated in Figs. EH and EHl and are: i) e~^e~ ZH, the Higgs strahlung or associated 
production with Z gauge bosons, ii) e+e" —>• Huu, the H^+H^” fusion production, and 
—>• He~^e~, the ZZ fusion production. Hi) e'^e~ —> tiH, the associated production 
with a tt pair. In addition, in the MSSM we also have: iv) e'^e~ —>■ h^A^,H^A^ and 
e’''e“ —> the pair production of two Higgs bosons, either neutral or charged, as 

illustrated in Fig. EH The cross sections of these processes as functions of the corresponding 
Higgs masses are illustrated in Figs. EH and EH for various center of mass energies (a/s). 
e’''e“ —> ZH and e^e“ —> Hui? are the leading production modes for a SM Higgs bosons. 
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FIG. 36; Higgs boson production in association with tt pairs at e+e colliders. 





FIG. 39: MSSM Higgs boson production cross sections in e'’"e collisions, for two choices of tan /3: 
tan/3 = 3 (solid) tan/3 = 30 (dotted), at -v/s = 350 GeV, for and H^, and at y^ = 800 GeV, for 


and H^. From Ref 


. Q. 
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Their relative size varies with the center of mass energy, since a{e'^e~ —>■ ZH) scales as 1/s 
(s-channel process), while a{e'^e~ —>■ Huu) scales as log(s) (t-channel process). —>■ Huu 

always dominates over e’''e“ —> He^e~ by almost one order of magnitude. In the MSSM, 
—> Huu plays a lesser role, due to the suppression of the VVH coupling {V = 1T^, Z), 
but the and H^H~ pair production modes become important. ttH production 

is always very rare, in particular at center of mass energies around 500 GeV or lower, but it 
plays a really important role at higher energies, around 800 GeV-1 TeV, for the determination 
of the top-quark Yukawa coupling, as we will discuss later. Other rare production modes 
that could play an important role in determining some Higgs boson properties are the double 
Higgs boson production modes: HHZ and HHuu 

With a LG running at energies between 350 GeV and 1 TeV, one or more Higgs bosons can 
be observed over the entire mass spectrum and all its properties can be precisely studied. 
Reconstructing the recoiling Hl~ mass (for I = e,fi) in HZ HHl~, where 

the Z is mono-energetic, allows an excellent and model independent determination of the 
Higgs boson mass. Z ^ qq decays can also be used and actually provide a very large 
statistics. Accuracies of the order of 50-80 MeV can be obtained, depending on the center of 
mass energy and the Higgs boson mass. The spin and parity of the Higgs boson candidate 
(expected to be = O’*") can be determined in several ways, among others: i) from the 
onset of a{e~^e~ —> ZH), since the energy dependence near threshold strongly depend on 
the quantum number of the radiated H (see Fig. EDI); a) from the angular distribution 
of H and Z in e’''e“ ^ ZH 4/; Hi) from the differential cross section in e+e“ —>■ tiH. 

Finally, a high energy LG will measure the Higgs boson couplings to unprecedented 
precision and in a model independent way. Thanks to the precise knowledge of the initial 
state energy conhguration (once the initial state radiation, or beam strahlung, has been 
properly taken into account), one can indeed measure both a{e^e~ HZ HHl~) and 
(T(e+e“ — W*W*h>u HuD), reconstructing the mass recoiling against the Hl~ or vu 
pair, and from there determine in a model independent way the isolated HZZ and HWW 
couplings^. This is probably the most important intrinsic difference between measuring the 
Higgs couplings at a lepton versus a hadron collider. At a lepton collider the Higgs couplings 

® We notice that the two production modes e+e” ^ HZ Hvv and e+e" ^ W*W*vv Hvv can be 
well isolated, since their distribution in the invariant vi> mass are very distinctive. See Ref. [l^ for further 
details. 
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FIG. 40: The e+e ZH cross 
spin = 0"'“, 1 ~, 2'’'. From Ref. 


section energy dependence near threshold for Mh 

□ 


120 GeV and 



FIG. 41: The theoretical predictions (bands, due to the uncertainty in the quark masses and q;<j) and 
experimental accuracies (dots with error bars) for the SM Higgs branching ratios at y/s = 350 GeV 
with 500 fb“^ of data. From Ref. 


to the weak gauge bosons, i.e. the Higgs couplings associated to the production mode (j/p of 
Section Fill D 2fl . can be isolated in a model independent way. Any other coupling can then 
be also determined in a model independent way, measuring the individual Br{H —> XX) 
in e’''e“ ^ HZ followed by if —> XX. Several recent studies have confirmed the possibility 
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FIG. 42: Higgs couplings determination at a LC (TESLA in this case) with 500 fb“^ of data, for 
= 120 GeV. The labels with fusion and w/o fusion refers to ITlF-fusion inclusion/exclusion. 
From Ref. 


u Cie V. 

f. 3- 


of determining Higgs couplings to both gauge boson and fermions within a few percent (2- 
5%). For instance, for a Higgs boson of = 120 GeV, the bottom-quark Yukawa coupling, 
yb, could be determined within 2%, the r one, within 5%, and the charm quark one, yc 
within 6% due to the larger error on the charm quark mass. This will test in a very stringent 
way the proportionality of the Higgs couplings to the mass of the interacting fermion. Even 
the indirect coupling of the Higgs boson to a pair of gluons, arising at the one-loop level 
(see Section. mH), will be determined with a precision of 4-5%. This will allow an indirect 
check of the top-quark Yukawa coupling for a SM Higgs boson (when the top-quark loop 
dominates), or will result in some anomalous coupling if new physics contributes in the 
loop. In this second case, probably, other Higgs coupling will show anomalous behaviors. 
The attainable precisions at a LC running at a/s = 350 GeV, and with a 500 fb“^ of data, 
are summarized in Fig. 

A LC will then be the ideal machine to discover small new physics effects. For instance, 
it could play a fundamental role in distinguishing the Higgs boson of the MSSM from 
the SM Higgs boson, in regions of the MSSM parameter space close to the decoupling limit. 
With this respect. Fig. 1^ shows two combined hts: the hrst one combines the ratio between 
the bottom-quark Yukawa coupling to and its SM expected value and the ratio between 








4 



Vs (GeV) 


FIG. 43: Cross section for e'’'e — > ttH, in the SM, as a function of the center of mass energy ^/s, 
for two different values of the Higgs boson mass. 

the coupling to W bosons and its SM expected value; the second one combines the same 
ratios for the bottom-quark and r-lepton Yukawa couplings. Different contours illustrates 
the precision with which the correlation between the combined ratios of couplings can be 
measured. On the same plots we see what the ratios would look like in the MSSM for 
different ranges of Ma- A definite distinction between SM and MSSM Higgs bosons can 
clearly be established. 

The only problem in completing a full study of the Higgs boson couplings is the deter¬ 
mination of the top-quark Yukawa coupling, and of the Higgs self couplings. The top-quark 
Yukawa coupling is indirectly determined by measuring the tiH cross section, when the 
Z contribution is under control. This cross section is very small at ^/s = 500 GeV, and 
peaks around y/s = 800 GeV (see Fig. m Recent studies show that with a LG operating 
at a/s = 500 GeV the top-quark Yukawa coupling, yt, for a Mh = 120 GeV Higgs boson, 
will probably be determined only at the 25% precision level, while with a LG operating at 
\/s = 800 GeV precisions as high as 5 — 6% becomes available. An updated summary of the 
existing studies can be seen in Fig. El The initial phase of a high energy LG will therefore 
not be able to give us probably the most important Yukawa coupling, and with this respect 
the role of the LHG becomes crucial, since, as we have seen in Section [111 D 21 the LHG can 
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H ^ bb semileptonic 
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H —> bb hadronic 


H ^ WW 2 like sign leptons 


H —> WW single lepton 


4 channels combined 



FIG. 44: Expected accuracies for the measurement of the top-quark Yukawa coupling, yt, in the 
process e+e" — > tlH, as a function of Mh, for y^=800 GeV and 1 ab“^ of data in various decay 
channels. From Ref. 

obtain yt within 10-15% accnracy, althongh with some intrinsic model dependence. 

The prodnction of Higgs boson pairs is also very rare (see Fig. 03]), and the measurement 
of the Higgs boson self-couplings will probably have to wait for a very high energy Linear 
Collider, like the CLIC collider, a multi-TeV machine being studied at CERN. 

Finally, the total width of a SM-like Higgs boson can be determined in a model inde¬ 
pendent way by using any well measured branching ratio. For example, one can use that 
T = T{H WW*)/Br{H —>• WW*), where Br{H —>• WW*) is measured directly and 

r(ff —>■ WW*) can be calculated from the direct determination of the HWW coupling. 

IV. HIGHLIGHTS OF THEORETICAL CALCULATIONS IN HIGGS BOSON 
PHYSICS 

I would like to conclude these lectures by reviewing some important theoretical results 
that have recently been obtained in the calculation of Higgs boson physics observables, i.e. 
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0.3 


SM Double Higgs-strahlung: e"^ e' ^ ZHH 



FIG. 45: Cross section for e"’“e“ ^ HHZ, in the SM, as a function of Mh, for two values of the 
center of mass energy ^/s. The dots with error bars show the achievable experimental accuracies 
for 1000 fb~^ of data. From Ref. [l^ . 

total and differential cross sections. Having to limit my discnssion to a selection of topics, I 
prefer to focns on hadron colliders, since the discovery of a Higgs particle very mnch depends 
on onr ability to provide precise theoretical predictions for hadronic cross sections. In this 
context I will concentrate on a few processes for which some ontstanding progress has been 
made recently. 

The cross section for pp and pp collisions to prodnce a hnal state containing a Higgs 
bosons {H) can be schematically written as: 

a{pp,pp ^ H + X) = Y^ J dxidx2fi{xi)ff’^{x2)a{ij ^ H + X) , (118) 

ij 

where the partonic cross section d{ij —>■ H + X) is convolnted with the Parton Distribntion 
Fnnctions (PDF) of partons i and j. ff’^ix) denotes indeed the PDF of parton i into a proton 
(or anti-proton) and can be interpreted as the probability of hnding parton i into a proton (or 
anti-proton) with a fraction x of its longitndinal momentnm. Both the partonic cross section 
and the parton distribntion fnnctions are calcnlated pertnrbatively. At hadron colliders, the 
most important effects arise from strong interactions, and it is therefore mandatory to have 
the QCD pertnrbative expansion of a{pp,pp ^ H + X) nnder control. At each order in 
the pertnrbative expansion, the calculation of both d^ij ^ H + X) and ff’^^x) contains 
ultraviolet divergences that are subtracted through a standard renormalization procedure. 
This, at each hnite order, leaves a dependence on the renormalization scale, pr. In the same 
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way, when the PDF’s are dehned, a factorization scale /ip is introduced in the calculation 
of f[’^{x). The dependence on both /i^ and /i^ is indicative of the residual theoretical 
uncertainty present at a given perturbative order, and should improve the higher the order 
of QCD corrections that are taken into account. 

Indeed, it is well known that the theoretical predictions for most Higgs production 
hadronic cross sections at lowest or leading order (LO) are affected by a very large renor¬ 
malization and factorization scale dependence. In general, at least the next-to-leading order 
(NLO) of corrections need to be calculated and this should stabilize or improve the theoret¬ 
ical prediction for the cross section, making the residual theoretical uncertainty comparable 
or smaller than the corresponding experimental precision. In some cases, as we will see, 
even next-to-next-to-leading order (NNLO) QCD corrections are necessary to obtain reli¬ 
able theoretical predictions. 

To be more specihc, the NLO cross section for pp,pp H + X can in full generality be 
written as: 


where we have made explicit the dependence on both renormalization and factorization scale, 
denote the NLO PDF’s, while is the parton level cross section calculated at NLO 


as: 


-NLO -LO I f--NLO 


where 6aij represents the 0{as) real and virtual corrections: 

^-NLO _ -ij I -ij 

^virt T ^real ' 


( 120 ) 


( 121 ) 


In a similar way, and with due differences, we could write A lot of theoretical effort 

has gone in recent years in the calculation of and for several Higgs production 

processes. I have collected in Table |1] all the existing work on higher order QCD corrections 
to Higgs production modes. My apologies for any omission! The result of this effort can 
be naively summarized by investigating the residual renormalization (pr) and factorization 
(Pf) scale dependence in each SM Higgs boson production mode. This is illustrated in 
Fig. uni for the case of the LHC. For each production mode the scale p = pr^= pp has been 
varied according to the corresponding original literature as indicated in Fig. (321 such that 
the scale interval can be different case by case. Since different production modes have very 
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TABLE I: Existing QCD corrections for various SM Higgs production processes. 


process 

0'NLO,NNLO by 

gg^H 

S.Dawson, NPB 359 (1991), A.Djouadi, M.Spira, P.Zerwas, PLB 264 (1991) 

C. J.Glosser et aL, JHEP 0212 (2002); V.Ravindran et aL, NPB 634 (2002) 

D. de Florian et al.^ PRL 82 (1999) 

R. Harlander, W.Kilgore, PRL 88 (2002) (NNLO) 

C.Anastasiou, K.Melnikov, NPB 646 (2002) (NNLO) 

V.Ravindran et aL, NPB 665 (2003) (NNLO) 

S. Catani et al. JHEP 0307 (2003) (NNLL), G.Bozzi et aL, PLB 564 (2003) hep-ph/0508068 

qq (IF, Z)H 

T.Han, S.Willenbrock, PLB 273 (1991) 

O.Brien, A.Djouadi, R.Harlander, PLB 579 (2004) (NNLO) 

qq qqH 

T.Han, G.Valencia, S.Willenbrock, PRL 69 (1992) 

T.Figy, C.Oleari, D.Zeppenfeld, PRD 68 (2003) 

qq, gg ttH 

W.Beenakker et al., PRL 87 (2001), NPB 653 (2003) 

S.Dawson et aL, PRL 87 (2001), PRD 65 (2002), PRD 67,68 (2003) 

qq, gg blH 

S.Dittmaier, M.Kramer, M.Spira, PRD 70 (2004) 

S.Dawson et aL, PRD 69 (2004), PRL 94 (2005) 

gb{b) ^ bib)H 

J.Campbell et aL, PRD 67 (2003) 

bb^ H 

D.A.Dicus et al. PRD 59 (1999); C.Balasz et aL, PRD 60 (1999). 

R.Harlander, W.Kilgore, PRD 68 (2003) (NNLO) 


different cross sections, we have separated them into two plots, containing the leading and 
sub-leading production modes respectively. The bands in Fig. EHlare in no way indicative of 
the overall theoretical error, since they do not include systematic errors coming from PDFs 
and other input parameters. Moreover the effect of setting 7 ^ /tf needs and has been 
investigated case by case, but it is not included in Fig. HHl Nevertheless Fig. HHl gives us a 
qualitative idea of the perturbative stability of the existing theoretical predictions for the 
SM Higgs boson production cross sections. Overall the existing theoretical predictions are 
in good control. 

n 

In a parallel series of lectures given at this school you have been exposed to the 
complexity of higher order QCD calculations, and to the variety of techniques that have 
been developed to perform them. I will not then directly proceed and comment about the 
results of some higher order calculations in Higgs physics. In particular, I would like to report 
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LHC, Vs = 14TeV, Mft/2 < n < 2M,, LHC, = 14TeV, /io = + Mft/2 




FIG. 46: Residual renormalization and factorization scale dependence {fJ^ = Hf) of the SM 

Higgs boson production cross section, when all available orders of QCD corrections are included. 
The scale /r is varied according to the original work present in the literature, and can therefore be 
slightly different case by case. 



FIG. 47: The gg ^ H production process at lowest order. 

about: i) the calculation of gg ^ H at NNLO of QCD, a pioneer effort that has provided 
for the hrst time a reliable theoretical results for the most important Higgs production 
mode at hadron colliders; ii) the calculation of pp,pp ttH,bbH at NLO of QCD, a 
challenging task, due to the many massive degrees of freedom involved, that has provided for 
the hrst time reliable theoretical predictions for the cross sections of these two physically very 
important production modes. Both i) and ii) rely on the development of several innovative 
techniques that have allowed the successful completion of both calculations. Given the 
degree of technicalities involved, I will not review them in detail, but only point to the kind 
of difficulties that had to be faced. The interested reader can hnd all necessary technical 
details in the original literature, listed in Table Q] and in the bibliography. 
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FIG. 48: The top-quark loop contribution to gg ^ H gives origin to a ggH effective vertex in the 
nit ^ oo limit. 


A. 


gg 


H at NNLO 


Most of the basic ideas that motivate the techniques used in the NNLO calculation 
of the cross section for the gg —y H production process have been already introduced in 
Section 1111 A 41 where we discussed the H ^ gg loop-induced decay. In particular we know 
that in the SM, the main contribution to gg H comes form the top-quark loop (see 
Fig. |17|) since: 


o'lo — 




( 122 ) 


288v/27r 

where Tq = and A^{Tq) < 1 with A^{Tq) —> 1 for Tq —>• oo. 

As we saw in Section Qll A 41 one can work in the inhnite top-quark mass limit and reduce 
the one-loop Hgg vertex to a tree level effective vertex, derived from an effective Lagrangian 
of the form: 

C,„ = Lc-(a,)G“'‘>'G;„ , (123) 

where the coefficient C{as), including NLO and NNLO QCD corrections, can be written as: 


C{,as) = - — 
3 TT 


1 I fOls 

1 + Cl-h C2 I - 

TT \ TT 


+ 


(124) 


NLO and NNLO QCD corrections to gg ^ H can then be calculated as corrections to the 
effective Hgg vertex, and the complexity of the calculation is reduced by one order of loops. 

The NLO order of QCD corrections has actually been calculated both with and without 
taking the infinite top-quark mass limit. The comparison between the exact and approximate 
calculation shows an impressive agreement at the level of the total cross section, and, in 
particular, at the level of the iC-factor, i.e. the ratio between NLO and LO total cross 
sections {K = o'nlo/ci'lo): as illustrated in Fig. HHl If is indeed expected that methods 
like the infinite top quark mass limit may not reproduce the correct kinematic distributions 
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FIG. 49: The NLO cross section for gg ^ H as a function of Mh- The two curves represent the 
results of the exact calculation (solid) and of the infinite top-quark mass limit calculation (dashed), 
where the NLO cross section has been obtained as the product of the iL-factor {K = (Tnlo/^lo) 
calculated in the rrit^oo limit times the LO cross section. From Ref. [l^ . 

of a given process at higher order in QCD, bnt are very reliable at the level of the total 
cross section, in particnlar when the cross section receives large momentnm independent 
contribution at the hrst order of QCD corrections. As for the H ^ gg decay process, the 
NLO corrections to gg ^ H are very large, changing the LO cross section by more than 50%. 
Since the gg ^ H is the leading Higgs boson production mode at hadron colliders, it has 
been clear for quite a while that a NNLO calculation was needed in order to understand the 
behavior of the perturbatively calculated cross section, and if possible, in order to stabilize 
its theoretical prediction. 

Recently the NNLO corrections to the total cross section have been calculated using the 
inhnite top-quark mass limit (see Table H]). The calculation of the NNLO QCD corrections 
involves then 2-loop diagrams like the ones shown in Fig. EOl instead of the original 3-loop 
diagrams (a quite formidable task!). Moreover, thanks to the 2 —1 kinematic of the 
gg ^ H process, the cross section has in one case be calculate in the so called soft limit, 
i.e. as an expansion in the parameter x = M‘jj/s about x=l, where s is the partonic center 
of mass energy (see paper by Harlander and Kilgore in Table Q}. The n-th term in the 
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FIG. 50: Two-loop diagrams that enter the NNLO QCD corrections to gg H. 



FIG. 51: iti-factor for gg ^ H at the LHC (\/s = 14 TeV), calculated adding progressively more 
terms in the expansion of Eq. (|120l) . From Harlander and Kilgore as given in Table HI 


expansion of of Eq. (HHl); 


CTi 


= E ? 


n>0 


cr, 


(n) 
ij ) 


( 126 ) 


can then be written in the soft limit (x —> 1) as follows: 


cr, 


in) 


2n-l 


= - x) + ^ 6, 


(n) 

k 


''ln^(l — x) 


k=0 


X 


00 2n—1 




( 126 ) 


J + 1=0 k=0 


purely soft terms 


collinear+hard terms 


where we have made explicit the origin of different terms in the expansion. The NNLO 
cross section is then obtained by calculating the coefficients , and cj^\ for / > 0 

and fc = 0,..., 3. In Fig. |^we see the convergence behavior of the expansion in Eq. (I12(i|l . 
Just adding the hrst few terms provides a remarkably stable iT-factor. The results shown in 
Fig. 1^ have been indeed conhrmed by a full calculation where no soft approximation 
has been used. 


The results of the NNLO calculation 




M Ml are illustrated in Figs. EH and 


In 
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^7= 14TeV 



FIG. 52: Cross section for gg ^ H at the LHC (\/i = 14 TeV), calculated at LO, NLO and NNLO 
of QCD corrections, as a function of Mh-, for gp = = Mh/4:. From Harlander and Kilgore in 

Table D 


a(pp ^H+X) [pb] V^= 2 TeV 


a(pp^H+X) [pb] V^= 14 TeV 



Mjj [GeV] 


FIG. 53: 
cross section for 
function of Mp. 
value gp 



Mjj [GeV] 


Residual renormalization/factorization scale dependence of the LO, NLO, and NNLO 
gg —> H, at the Tevatron (-y/i = 2 TeV) and at the LHC {^/s = 14 TeV), as a 
The bands are obtained by varying gpi^ = gp hy a factor of 2 about the central 
gR = MH/4:. From Ref. 


U Ulcll lT 

f. Q 


Fig. ESI we can observe the convergence of the perturbative calculation of u{gg H), since 
the difference between NLO and NNLO is much smaller than the original difference between 
LO and NLO. This is further conhrmed in Fig.Eni where we see that the uncertainty band of 
the NNLO cross section overlaps with the corresponding NLO band. Therefore the NNLO 
term in the perturbative expansion only modify the NLO cross section within its NLO 
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FIG. 54; NNLL and NNLO cross sections for Higgs boson production via gluon-gluon fusion at 
both the Tevatron and the LHC. From Ref. Q 


theoretical uncertainty. This is precisely what one would expect from a good convergence 
behavior. Moreover, the narrower NNLO bands in Fig. 1^21 shows that the NNLO result is 
pretty stable with respect to the variation of both renormalization and factorization scales. 
This has actually been checked thoroughly in the original papers, by varying both fiji and 
Pi? independently over a range broader than the one used in Fig. EHl 

The NNLO cross section for gg ^ H has been further improved by Catani et ah Q by 
resumming up to the next-to-next-to leading order of soft logarithms. Using the techniques 
explained in their papers, they have been able to obtain the theoretical results shown in 
Figs. EH and EH for the total and differential cross sections respectively. In particular, we 
see from Fig. I^that the NNLO and NNLL results nicely overlap within their uncertainty 
bands, obtained from the residual renormalization and factorization scale dependence. The 
residual theoretical uncertainty of the NNLO+NNLL results has been estimated to be 10% 
from perturbative origin plus 10% from the use of NLO PDF’s instead of NNLO PDF’s. 
Moreover, in Fig. EH we see how the resummation of NNL crucially modify the shape of the 
Higgs boson transverse momentum distribution at low transverse momentum (qt), where the 
soft \n{Mjj/q^) are large and change the behavior of the perturbative expansion in ag 


B. pp,pp —> ttH and pp,pp bbH at NLO 

The associated production of a Higgs boson with heavy quark pairs, pp,pp —>■ tiH and 
bbH, has been for a while the only Higgs production process for which the NLO of QCD 
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FIG. 55: The qt spectrum at the LHC with Mh = 125 GeV. The upper plots show: (left) setting 
= fJ-F = Q = Mh, the results at NNL+LO accuracy compared with the LO spectrum and the 
finite component of the LO spectrum; (right) the uncertainty band from the variation of the scales 
fiji and lip Sit NLL+LO accuracy. The lower plots show the same at NNLL+NLO accuracy. From 


Ref. 


56|. 


corrections had not been calcnlated. Given the relevance of both prodnction modes to Higgs 
physics (see discnssion in Section ITTTdI) and the large renormalization and factorization 
scale dependence of the LO cross sections (see, e.g., the LO cnrves in Figs. IHHl and . a 
full NLO calculation was mandatory. This has been completed in the papers by Beenakker 


et ah j^, 1^ and Dawson et ah 


59 


61 


62| | listed in Table HI and we will briefly report 


about their most important results in this Section. 
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FIG. 56: Tree level Feynman diagrams for qq —> ttH and gg —> ttH. Analogous diagrams, with 
t ^ b, contribute to the qq —> hbH and gg —> bbH processes. 


The NLO calculation of pp,pp QQH (for Q = b,t) presents several challenges, since it 
has to deal with a 2 —>■ 3 process that involves all massive particles in the final state. At 
tree level, QQH proceeds through the qq, gg —> QQH parton level processes illustrated in 
Fig. EElfor Q = t. As expected, the qq QQH dominates for large fraction of the parton 
longitudinal momentum x, while gg QQH dominates at small x. This translate into the 
fact that the parton level cross section for tiH production is dominated by qq tiH at 
the Tevatron, and by gg tiH at the LHC, while bbH production is always dominated 
by gg —^ bbH. The 0{as) virtual corrections include up to pentagon diagrams, such that 
the problem of calculating scalar and tensor integrals with up to five denominators, several 
of which massive, has to be faced. Most integrals have to be calculated analytically since 
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FIG. 57: Example of ©(a^) virtual corrections to gg —> tiH: pentagon diagrams. The circled 
crosses denote all possible insertion of the final Higgs boson leg. All t-channel diagrams have 
corresponding u channel diagrams, where the two initial state gluon legs are crossed. Analogous 
diagrams with t ^ b contribute to the 0{as) virtual corrections to gg —> bbH. 



FIG. 58: Examples of 0(as) real corrections to gg tiH. The circled crosses denote all possible 
insertion of the final Higgs boson leg. Analogous diagrams with t —> 6 contribute to the 0{as) real 
corrections to gg —> bbH. 


both ultraviolet (UV) and infrared (IR) singularities have to be extracted. The 0{as) real 
corrections involve factoring out IR (soft and collinear) divergences from a 2 —>■ 4 phase 
space with several massive particles. Samples of Feynman diagrams corresponding to the 
0{as) virtual and real corrections are illustrated in Figs. EZI and EH Several new methods 
and algorithms have been used by the two collaborations that have calculated the NLO 
cross section for pp, pp —> QQH, and we refer to their papers for all technical details. In the 
following we will comment separately on tiH and bbH NLO results. 

For tiH production, the most important outcome of the full NLO calculation is illus¬ 
trated in Fig. EH where the renormalization (pr) and factorization scale (pf) dependence 
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FIG. 59; Dependence of (Tlo,nlo{pP-,PP ttH) on the renormalization/factorization scale fi, 
at both the Tevatron and the LHC, for Mh = 120 GeV. The reference scale /Un is taken to be 
^Q = mt + Mh/2. From Refs. 


6C| and 


y. 


of the LO and NLO total inclusive cross section is presented, for a SM Higgs boson mass 
of Mh = 120 GeV, at both the Tevatron and the LHC. We note that the factorization and 
renormalization scales have been set equal in the plots of Fig. EHl Rr = Rf = R, while in the 
original work both scales have been hrst varied independently to verify that /xr = /if = R is 
not a particular point at which both scale dependences accidentally mutually cancel. It is 
evident from the plots of Fig. EHl that the NLO total cross section sensitivity to /xr and pp 
is drastically reduced with respect to the LO cross section. Indeed, the residual systematic 
error coming from scale dependence is at NLO reduced to about 15%, as opposed to more 
than 100% of the LO cross section. The NLO predictions for ttH production can now be 
conhdently used to interface with experimental analyses. 

Let us now consider bbH production®. Naively, one would expect the calculation of bbH 
production at NLO to follow that of tiH, with the universal replacement of the top-quark 
mass with the bottom-quark mass, mb. However, the theoretical prediction of bbH 

production at hadron colliders involves a few subtle issues not encountered in the calculation 
of tiH production. Indeed, both from the experimental and theoretical standpoint, it is 
important to distinguish between inclusive and exclusive bbH production. More specifically, 
the production of a Higgs boson with a pair of b quarks can be detected via: (i) a fully 

® For an updated review see Ref. 0 
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exclusive measurement, where both h jets are observed, (u) a fully inclusive measurement, 
where no b jet is observed, or [Hi) a semi-inclusive measurement, where at least one b jet is 
observed. 


Experimentally, b quarks are identihed or tagged by imposing selection cuts on their 
transverse momentum and their angular direction with respect to the beam axis or pseudo¬ 
rapidity. Inclusive modes have larger cross sections, but also larger background, such that 
more exclusive modes are often preferred experimentally. Moreover, only the exclusive and 
semi-inclusive modes are unambiguously proportional to the bottom-quark Yukawa coupling. 

Theoretically, different approaches may be adopted depending on the fact that a hnal state 
b quark is either treated inclusively (untagged) or exclusively (tagged). Indeed, when a hnal 
state b quark is not identihed through some selection cuts, the corresponding integration 
over its phase space, in particular over its transverse momentum, gives rise to logarithms of 
the form: 



( 127 ) 


where and fin represent the lower and upper bounds of the integration over the transverse 


momentum of the hnal state b quark, fin is typically of 0{Mh) and therefore, due to the 
smallness of the bottom-quark mass, these logarithms can be quite large. Additionally, the 
same logarithms appear at every order in the perturbative expansion of the cross section in 
as, due to recursive gluon emission from internal bottom-quark lines. If the logarithms are 
large, the convergence of the perturbative expansion of the cross section could be severely 
hindered and it can be advisable to reorganize the expansion in powers of q;"A™, further 
resumming various orders of logarithms via renormalization group techniques^. 

Currently, there are two approaches to calculating the inclusive and semi-inclusive cross 
sections for Higgs production with bottom quarks. Working under certain kinematic ap¬ 
proximations, and adopting the so-called five-flavor-number scheme (SENS), the collinear 
logarithms, A^, can be factored out and resummed by introducing a bottom-quark Parton 
Distribution Function (PDF) |^Q,|^. This approach restructures the calculation to be an 
expansion in both as and A^^. At tree level, the semi-inclusive production is then described 
by the process bg bH illustrated in Fig. 1^01 while the fully inclusive production process 

^ The logarithms mentioned here also appear in the ttH calculation but, since gu is typically of the order 
of TTit, the logarithms are small and the convergence of the perturbative expansion in as is preserved. 


84 



becomes bb H, illustrated in Fig. Alternatively, working with no kinematic approx¬ 
imations, and adopting the so-called four-flavor-number scheme (4FNS), one can compute 
the cross section for pp, pp bhH at hxed order in QCD with no special treatment of the 
collinear logarithms, considering just the parton level processes qq,gg —> bbH illustrated in 
Fig. inni (with t ^ b) and their radiative corrections. 



FIG. 60: Tree level Feynman diagram for bg bH in the 5FNS. 



h 


FIG. 61: Tree level Feynman diagram for bb ^ H in the 5FNS. 

The fully exclusive bbH production cross section can only be computed in the 4FNS 
framework. As far as the inclusive and semi-inclusive production cross sections go, the 
comparison between the 4FNS and 5FNS needs to consider QCD corrections beyond the 
LO, in order to work with stable results. Indeed, the two calculation schemes represent 
different perturbative expansion of the same physical cross section, and therefore should 
agree at sufficiently high order. The discussion to follow is based on the NLO calcu 
qq^ gg ^ bbH 0| and bg —> bH Q, and on the NNLO calculation of bb ^ H 

It should be noted that our discussion for the production of a scalar Higgs boson with 
bottom quarks applies equally well to the production of a pseudoscalar Higgs boson. In 
fact, if one neglects the bottom-quark mass in the calculation of the NLO corrections, the 
predictions for bhA^ is identical to those for bbh^{H^) upon rescaling of the Yukawa couplings 
(see Section FH D 4fl . On the other hand, for massive b quarks, the situation becomes more 
complicated due to the 75 matrix appearing in the bbA^ Yukawa coupling. The 75 Dirac 
matrix is intrinsically a four-dimensional object and care must be taken in its treatment 


ation of 
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when regularizing the calculation in dimensional regularization (d 7 ^ 4). However, bottom- 

2 _ 

quark mass effects are expected to be small, 0{^), and predictions for bbh^, upon rescaling 
of the Yukawa coupling, provide good estimates of bbA^ production even in the massive 
6 -quark case. In the following we will present results for the exclusive, semi-inclusive, and 
inclusive cross sections separately. 

The fully exclusive bbH NLO total cross section is illustrated in Figs. IHH and 1221 both 
for the Tevatron and for the LHC. Both b quarks in the final state are identified by the 
cuts explicitly given in Fig. |22l which have been chosen to closely mimic experimental 
searches. The curves in Fig. 1221 show the dependence of the LO and NLO exclusive bbH 
cross section from both renormalization and factorization scales (set equal in these plots, i.e. 
IJ'R = IJ'F = yb)- The two sets of curves represent the case in which the bottom-quark mass in 
the bottom-quark Yukawa coupling is renormalized in the on-shell scheme (OS', blue curves) 
or in the modified Minimal Subtraction scheme (MS', red curves). As we have already 
observed in Section UlI A 31 when we considered the decays of a Higgs boson into quark pairs, 
the bottom-quark renormalized mass varies substantially when the renormalization scale /i 
is varied from scales of the order of mb to scales of the order of Mh, and therefore it is 
important to know how the large logarithms that determine the running of are treated 
in the different renormalization schemes. This is particularly true for the factor mb that 
appears in the bottom-quark Yukawa coupling, yb = mb/v, since the cross section depends 
quadratically on yb- It is much less relevant for the cross section kinematic dependence on 
mb, coming from the amplitude square or from the integration over the final state phase 
space. For this reason, the OS and MS labels in Fig. |221 refer to the cases in which only 
the bottom-quark mass in yb is renormalized one way or the other. The kinematic 6 -quark 
mass is always taken to be the pole mass. Looking at the curves in Fig. |221 we learn that, 
as expected, the NLO set of QCD corrections stabilizes the total cross section, drastically 
reducing the scale dependence of the LO cross section. Moreover we see that there is a non 
negligible dependence on the bottom-quark mass renormalization scheme. This dependence 
is intrinsic of any perturbative calculation, and should decrease the more orders are added in 
the perturbative expansion of a given physical observable. However, from the behavior of the 
residual scale dependence, it is possible to estimate which renormalization scheme provides 
a better perturbative expansion. In the case illustrated in Fig. |221 both OS and MS NLO 
cross sections show a well defined plateau region where the cross section is very mildly 








FIG. 62: Scale dependence of the LO and NLO cross sections for pp — > bbH (Tevatron) and 
pp bbH (LHC), for Mh = 120 GeV. The curves labeled (Tlo,OS and (Tnlo,os use the OS 
renormalization scheme, while the curves labeled (Jlo,ms and (Tnlo,ms use the MS renormalization 
scheme. From Ref. y. 


dependent on the scale p. Nevertheless, the MS cross section overall performs better, since 
it has a more regular behavior also at small scales. As a result, MS is often the preferred 
choice in all processes that depend on the bottom-quark Yukawa coupling. Finally, Fig. IHHl 
shows the dependence of the total cross section from the Higgs boson mass, in both the SM 
(solid black curve) and the MSSM (dashed colored curves, corresponding to different values 
of tan/?), over a significant Mh range for both the Tevatron and the LHC. The quadratic 
growth of the MSSM cross sections with tan/? is evident, and this graphically confirms 
the possibility of finding evidence of new physics already at the Tevatron, if an MSSM-like 
2HDM with large tan/? is realized in nature (see Sec. 1111 iT^ . 

Let us now turn to the semi-inclusive and inclusive bbH cross sections. For both 
production modes, much effort has been spent recently in understanding the difference 
between 4FNS and 5FNS in order to assess the reliability of the existing theoretical 
resultsj^,0,13, Ell- First of all, let us briefly review the idea behind the introduction of a 
bottom-quark PDF, which naturally leads to the 5FNS framework. For the purpose of illus¬ 
tration, consider the prototype case depicted in Fig.lHH one of the hnal state bottom quarks 
is directly originating from the g ^ bh splitting of one of the initial state gluons, while the 
shaded blob represents all the possible non collinear conhgurations of the remaining parti- 
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FIG. 63: NLO cross section for pp bbH (Tevatron) and pp bbH (LHC) as a function of Mh, 
in the SM and in the MSSM with tan/3 = 10, 20 and 40. For the Tevatron the production process 
considered is pp —> bbhP with = 100,110,120, and 130 GeV, while for the LHC it is pp —> bbH^ 


with Mfjo = 120, 200,400,600, and 800 GeV. From Ref. 


6 £|. 



FIG. 64: Tree level Feynman diagrams for gg bbH illustrating the almost collinear emission of 
a bottom quark (upper leg). 


cles. In the ^ 0 limit, the g ^ bb conhguration gives origin to collinear singularities, 
when the two b quarks are emitted in the same direction of the splitting gluon. In the case of 
bbH production, the singularities appear in the —>• 0 phase space region, where is the 
transverse momentum of the upper leg bottom quark in Fig. If we take ruf, ^ 0, these 
singularities are regulated by the non zero 6-quark mass, leaving behind collinear logarithms 
Afo of the form given in Eq. (1I2II). The contribution to the total partonic cross section from 
this diagram can be written as: 


da. 


gg^bbH 


d^Pb 1 d?Pb 1 dr'pH 1 


-Y.\M 


gg^bbH I 


{2)2Ei2E2 (27r)3 2Eb (27r)3 2Eb (27r)3 2Eh 
■{2TiY5^'^\qi + q2-pb-Pb-PH) , 


2 


(128) 













where we have used the momentum notation of Fig. such that Ei ^2 are the energies of 
the initial gluons while E^i ^ are the energies of the hnal state particles. The amplitude for 
this process is denoted by M-gg^hln- Parameterizing the 6-quark propagator momentum as: 

= ZQi + f3n^ -I- , (129) 


where z is the fraction of the initial gluon momentum carried by the b quark, and fcj) is 
transverse to both gf and for an arbitrary vector n^, with k'^k±n = —k^, one can 
explicitly show (after some Dirac algebra manipulations) that: 

\-^gg^bbH\‘^ — ^ -Pq9(^) ^ ) ( 130 ) 

where Pqg{z) is the Altarelli-Parisi splitting function for g qq: 

P,, = + {1 - zf] , (131) 


Cp = (IV^ — l)/2/iV (for N = 3 colors), and M.gi^ih is the amplitude for gb bH. The 
previous equation is approximate, i.e. it has been obtained by neglecting higher powers of 
keeping only the terms that would give the most singular or leading contribution upon 
integration over the collinear 6-quark phase space: 

(Ppb 1 1 


dzdk"^ 


(132) 


(2ir)3 2£(, 16ir2 (1 - 2 ) ’ 

where we have used that in the small kp limit the transverse part of the four momentum of 
the outgoing quark, pj coincides with kp- Inserting Eqs. ()ld()|l and ()ld2|l into Eq. ()128|1 one 
hnds: 


dk'^ 


0.0 


d^gg^hbH — j^2 r)^^ggi,^^^^gb^bH 
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(133) 


The integration over k"^, with lower bound ml and upper bound /i^, gives origin to the 
collinear logarithm A;, introduced in Eq. (iniD. Moreover, when one convolutes with the 
gluon PDFs, g{xi,p) (for i = 1,2), of the two initial gluons to obtain the hadronic cross 
section, a bottom-quark PDF of the form: 

asipr) 


b{x,pF) = 


27r 


E y \yj 


(134) 


naturally appears. In Eq. (uni fip represents the factorization scale. The collinear Af, 
logarithms are factored out and then resummed in the bottom-quark PDF when the factor¬ 
ization scale is set to pF = yH- Indeed, Eq. (I134|l gives the bottom-quark parton density at 








the lowest order in as, while the leading logarithms are resummed via the DGLAP 

equation upon evolution: 


~ Jx y \yj 


TT 


(135) 


dlog/i 

The 5FNS approach is therefore based on the approximation that the outgoing b quarks 
are at small transverse momentum, since this is the region of phase space that is emphasized 
by the /ct expansion. The incoming b partons are given zero momentum at leading order, and 
acquire transverse momentum at higher order. With the use of a 6-quark PDF, the 5FNS 
effectively reorders the perturbative expansion to be one in and A^^. To see how this 
works, let us consider the perturbative expansion of the inclusive process 66 —> iP (Fig. IF)T|l 
which, according to what we just saw, is intrinsically of order a'^Af. At NLO, the virtual 
and real corrections to the tree level process make contributions of (^(afA^). However, at 
NLO, we must also consider the contribution from bg —> bH where the hnal state 6 is at 
high transverse momentum. This process makes a contribution of order a^Af, and is, thus, a 
correction of (P(A^^) to the tree level cross section. Similarly, at NNLO, besides the myriad 
of radiative corrections of (P(q;^A^), we must also include the contribution from the process 
gg —> bbH, where both 6 and 6 are at high pt- The contribution from these diagrams 
are of order a?, and are, thus, (P(A^^) (or NNLO) corrections to the tree level process 
66 —> H 1^ . The above discussion for bb ^ h also applies to the perturbative expansion 
of bg —>■ bH. In this case, the tree level process is of order a'^Aj, and the contribution from 


gg —bbH is a NLO correction of 0{A 


-l^ 


Q. 
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The comparison between 4FNS and 5FNS has been initially performed in the SM 
0, 0,[2ll. This has been crucial to understand several important issues. However, since 
the production of a Higgs boson with bottom quarks will only be physically interesting if 
the bottom-quark Yukawa coupling is enhanced beyond its SM value, we will present the 


Q. Figs. 


and ESI show the 


4FNS vs 5FNS comparison in the MSSM, with tan/3 = 40 
comparison for the inclusive and semi-inclusive total cross sections respectively, at both 
the Tevatron and the LHC. They both give the cross section as a function of the Higgs 
boson mass, chosen to be at the Tevatron and H^ at the LHC. In the inclusive case. 
Fig. ESI illustrates both the NLO and NNLO predictions. Both at LO and NLO (NNLO) the 
band are obtained by varying the renormalization and factorization scales independently as 
explained in the hgure captions. In both the inclusive and semi-inclusive case there is good 
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agreement between the 4FNS and 5FNS results within their respective scale uncertainties, 
although the 5FNS tends to always give slightly higher results. 
cnlo [ph] Tevatron, ^/s = 1.96 TeV unlo [ph] LHC, ^/s = 14 TeV 



quark jet identified in the final state in the 4FNS (at NLO) and 5FNS (at NLO and NNLO) as a 
function of the light and heavy MSSM Higgs boson masses, at both the Tevatron and the LHC. 
The error bands have been obtained by varying the renormalization (p,.) and factorization (pj) 
scales separately between po/4 and po (with ^o = mh + Mh/2) in the 4FNS, while keeping ^r = Mh 
and varying pj between O.IMh and 0.7Mh in the 5FNS (see Ref. for details). From Ref. fl¬ 


it is actually very satisfactory that the NLO (and NNLO) calculations for the semi- 
inclusive and inclusive bbH production agree within their systematic errors. Indeed, if all 
perturbative orders were to be considered, the two approaches would produce the very same 
outcome. The truncation of the perturbative series at a given order gives origin to discrep¬ 
ancies, because the 4FNS and 5FNS perturbative series are differently ordered. However, 
we see that considering the first (second) order of corrections already brings the agreement 
between the two schemes within the respective theoretical uncertainties. The results shown 
in Figs, inni and insi represent indeed a major advancement, since the comparisons existing 

671 1681 Rfll [7l| showed a pronounced disagreement. This was 


in the literature before Refs. 


mainly due to the absence of the (now available) NLO results for qq, gg —> bbH. 
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cnlo [pb] Tevatron, ^/s = 1.96 TeV 




100 105 no 115 120 125 130 150 200 250 300 350 400 450 500 

Mho [GeV] Mho [GeV] 

FIG. 66: Total NLO cross section in the MSSM for pp,pp —> b(b)h production at the Tevatron and 


the LHC as a function of We varied pr and pj independently to obtain the uncertainty 

bands, as explained in the text. The solid curves correspond to the 4FNS, the dashed curves to the 


5FNS. The error bands have been obtained by varying the renormalization and factorization 


{pf) scales separately between 0.2/ro and /tq (with po = mb + Mh/‘2)- From Ref. 
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